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1. Introduction 

^. 

O . Let a, b and c be positive integers and consider a hexagon with side lengths a,b,c,a,b,c 

whose angles are 120° (see Figure |l]). 

The subject of our interest is rhombus tilings of such a hexagon using rhombi with all sides 

' Y^ [ of length 1 and angles 60° and 120°. Figure |2| shows an example of a rhombus tiling of a 

^ ' hexagon with a = 3, 6 = 5 and c = 4. 

A first natural question to be asked is how many rhombus tiling of a fixed hexagon exist. 

A well known bijection between such rhombus tilings and plane partitions contained in an 

ax bx c box [P] and MacMahon's enumeration of plane partitions |11, Sec. 429, (7 — > 1; proof 



0^ 



X 



in Sec. 494] give the following answer: The number of all rhombus tilings of a hexagon with 



d ' side lengths a,b,c,a,b,c equals 



c 



i=ij=ik=i •' 

On a next level, one may ask for the number of rhombus tilings with special properties. In 
this paper, we address this question. We study rhombus tilings of a hexagon which contain 
the central rhombus and rhombus tilings of a hexagon which contain the 'almost central' 
rhombus above the centre. By the 'central rhombus' we mean the rhombus whose centre 
is equal to the centre of the hexagon (see Figure |3|, where the central rhombus is marked; 
furthermore, the example of a rhombus tiling in Figure |2| contains the central rhombus). By 
the 'almost central' rhombus above the centre we mean the horizontal rhombus whose lowest 
vertex is the centre of the hexagon (see Figure |^ . 

The main results of this paper are the following two theorems. 
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Theorem 1. Let a,b,c be positive integers with a = b {mod'2) and a ^ c {mod2). Then 
the number of rhombus tilings of a hexagon with side lengths a,b,c,a,h,c which contain the 
rhombus in the centre is 
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in case that a is odd, and 
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(1.3) 



m case that a is even, where the Pochhammer symbol (a)^ is defined by (a)^ := a(a+l) . . . (a+ 
k — 1) if k > 1 and (a)o := 1. 



The special case a = 6 in Theorem || was previously derived in |]|, Theorem 1,2]. (In order 
to see that the two sums in Theorem |l| are equal to two sums in Theorem 1 and Theorem 2 
from [y in this special case, one has to apply Bailey's transformation (|5.3| ) for balanced 
4F3-series.) 

The assumption about the parity of the side lengths a,b,c, i.e., that not all side lengths 
of the hexagon have the same parity, comes from the fact that this condition is necessary 
for the existence of a central rhombus in a rhombus tiling of a hexagon with side lengths 
a,b,c,a,b,c. There is also a result similar to Theorem || for a hexagon whose side lengths have 
the same parity. In this case we choose for the fixed rhombus the horizontal rhombus whose 
lowest vertex is the centre of the hexagon. (We could have also chosen the horizontal rhombus 
whose uppermost vertex is the centre of the hexagon). The number of these rhombus tilings 
is given by formulas quite similar to those in Theorem |l|, and, furthermore, the proofs of these 
formulas are analogous to the proofs of the formulas in Theorem IT]. We obtain the following: 



Theorem 2. Let a,b,c be positive integers with a = b = c {mod 2). Then the number of 
rhombus tilings of hexagon with side lengths a,b,c,a,b,c which contain the horizontal rhombus 



ENUMERATION OF RHOMBUS TILINGS 

whose lowest vertex is the centre of the hexagon is equal to 

[ " ^ "" i + j + k-l\ (l)c / &+C-2 X . g+b+e-l N 

l^nnn^T7Ti^j(6Ti)^l4Jl 4 r 
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in case that a is odd, and 
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in case that a is even. 



The special case a = 5 in Theorem |2| was previously derived in [Q, Theorem 1, 2]. (In 
order to see that the two sums in Theorem ^ are equal to the two sums in Theorem 1 and 
Theorem 2 from Q for this special case, one first has to apply a certain contiguous relation 
on the balanced 4-F3~series in Theorem ^. In each case we obtain the sum of a balanced 
3-F2-series, which can be evaluated by Saalschiitz's summation formula (see O, (2.3.1.3)]), 
and another balanced 4-F3-series. An application of Bailey's transformation to the latter 4-F3- 
series finally shows the equivalence of the result in Theorem and the result in P, Theorem 
1, 2] for the special case.) 

These enumeration results for rhombus tilings with a fixed rhombus are not only interesting 
because they add more results to the growing set of results on the enumeration of rhombus 
tilings with special properties. These enumeration also contribute to the interesting question 
of what a 'typical' rhombus tiling of a hexagon looks like. Cohen, Larsen und Propp address 
this question in P] . We are able to deduce the following two theorems from Theorem ffl and 
Theorem y. They confirm Conjecture 1 in Q for a special case. 

Theorem 3. Let a, (3 and 7 he non-negative real numbers. Then the probability that a 
rhombus tiling of the hexagon with side lengths a, b, c, a, b, c, where a ~ aN , b ~ (iN , c ~ 7A^ 
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and a = b^ c {mod2), contains the rhombus in the centre is asymptotically 




"^ ^ (1.6) 



as N tends to infinity. 

Theorem 4. Let a, (3 and 7 be non-negative real numbers. Then the probability that a rhom- 
bus tiling of the hexagon with side lengths a, 6, c, a, 6, c, where a ~ aN , b ~ PN, c ~ 7A'' and 
a = b = c {mod2), contains the 'almost central' rhombus above the centre is asymptotically 




"^ ^ (1.7) 



^^ (/3 + 7)(a + 7) 
as N tends to infinity. 

Roughly speaking, Conjecture 1 in [^ predicts the fohowing: Fix an arbitrary point (x, y) in 
the hexagon with side lengths a,b,c,a,b,c. Then the probability that a random rhombus tiling 
of a hexagon contains a vertical rhombus at this point tends to Va,/3,-y{x,y) as the hexagon 
becomes large. Here, 'Pa,/3,7 is a certain function (defined in |2|, Theorem 1]) that depends on 
the proportions a,/3,7 of the side lengths of the hexagon (see Theorem 0). Theorem H and 
Theorem § confirm this conjecture for the centre of the hexagon. (The reader should know 
that Cohen, Larsen and Propp use another coordinate system in |Q] than we do and therefore 
they are in the position to speak of vertical rhombi. But of course, their result can easily 
be translated into our coordinate system where we fix horizontal rhombi instead of vertical 
rhombi.) 

We want to direct the reader's attention to a side result of the present work, given in 
Lemma g. It expresses, for any fixed rhombus, the number of rhombus tilings containing that 
fixed rhombus as a triple sum. I am currently pursuing an asymptotic analysis of this triple 
sum with the ultimate goal of effectively proving Cohn, Larsen and Propp's conjecture. 

In order to prove Theorem ^ and Theorem ^ we make use of a bijection between rhombus 
tilings which contain the central rhombus, respectively rhombus tilings which contain the 
'almost central' rhombus above the centre, and non-intersecting lattice paths. This bijection 
actually works for arbitrary rhombi, i.e., the fixed rhombus need not be placed in the centre 
or next to the centre. The description of this bijection is the subject of Section ^. Since the 
number of non-intersecting lattice paths is given by a determinant due to Lindstrom, Gessel 
and Viennot (see Theorem ^), this bijection reduces the problem to evaluating a certain 
determinant with binomial entries, see Lemma [l|. In case that the fixed rhombus is situated 
in the centre or next to the centre we are able to evaluate the determinant, see Lemma |^ and 
Lemma H in Section B. We partly make us of a method that has already produced evaluations 
of other binomial determinants (see e.g. y], yj). In the course of evaluating the determinant 
we need an alternative expression for the number of rhombus tilings that we are interested 
in, in form of the aforementioned triple sum. This is given in Lemma § in Section y. Finally, 
in Section |5| we provide the proofs of Theorem |3| and Theorem Q. 

2. From rhombus tilings to non-intersecting lattice paths and determinants 

As mentioned before, the subject of this section is the bijection between non-intersecting 
lattice paths and rhombus tilings. In our context, non-intersecting lattice paths are disjoint 
paths in the lattice Z^ with steps in direction (1, 0) and (0, —1). Figure shows an example of 
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such a family of non-intersecting lattice paths. It consists of three paths, each one connecting 
an initial point Ai with a destination point Ei, i = 1,2, 3. 

Figure ^ and Figure |6| illustrate the bijection between rhombus tilings of a hexagon with 
side lengths a,b,c,a,b,c and non-intersecting lattice paths with initial points 

Ai = {i-l,c + i-l) (2.1) 

and destination points 

Ei = {b + i-l,i-l), (2.2) 

i = 1,2, ... ,a. Figure g shows the rhombus tiling from Figure |2| and an indication of the 
corresponding non-intersecting lattice paths. The initial points A'- and the destination points 
E'- of the paths in this figure are the centres of the sides of the rhombi that form the two 
sides of the hexagon with side lengths a. Roughly speaking these paths just describe 'the 
way down' from A'- to E'- on the three dimensional pile of cubes — elsewhere called plane 
partitions — which the rhombus tiling, when interpreted as three-dimensional object, gives. 

In order to obtain the corresponding family of non-intersecting lattice paths in Figure ^, 
we only have to change the 120° angles of the paths in Figure ^ into right angles. 

Therefore we already have a bijection between all rhombus tilings of a hexagon with side 
lengths a,b,c,a,b,c and all non-intersecting lattice paths starting in Ai and stopping in E'j. 
However, we are interested in special rhombus tilings — namely in those with one fixed 
rhombus. Because of that we have to look for a simple description of the non-intersecting 
lattice paths that correspond to the rhombus tilings with a fixed rhombus. Studying Figure S 
and Figure y again, we see that the marked rhombus in Figure ^ corresponds to the edge 
(x — l,y) — > {x,y) in Figure ^ (in our example x = 4 and y = 3). Therefore the rhombus 
tilings which contain a fixed horizontal rhombus correspond to the families of non-intersecting 
lattice paths which contain a fixed horizontal edge. Thus, we have found a bijection between 
rhombus tilings of a hexagon with side lengths a,b,c,a,b,c which contain a fixed horizontal 
rhombus and non-intersecting lattice paths with initial points Ai and destination points Ei, 
i = 1,2, . . . ,a, which contain a fixed horizontal edge. As we will see in a moment, the latter 
can be counted with the help of the following main theorem on non-intersecting lattice paths: 

Theorem 5 (Lindstrom-Gessel-Viennot). Let Ai,A2, . . . ,Ar and Ei,E2, . . . ,Er be points in 
T? . Then the determinant 

det {\P{A,^E,)\) (2.3) 

l<«,j<r 

is equal to 

5^sgna-|P+(A^K)|, 

where \P [Ai — > Ej)\ denotes the number of lattice paths connecting Ai to Ej, and {P'^ {A — > E^ 
denotes the number of families {Pi, P2, . . . , Pr) of non-intersecting lattice paths, the ith path 
Pi connecting Ai to E^ii\, 1 <i < r. 

This theorem can be found in ^, Theorem 1] or, alternatively, in jl^. Lemma 1]. 

Remark 1. Usually the theorem is applied to the following situation: There exists only 
one permutation a (normally the identity permutation) such that the set of families of non- 
intersecting lattice paths connecting Ai to J^g- (j), 1 < i < r, is not empty (e.g., this is the case 
for the points Ai and Ei defined in (^]^) and (2.2)). Then Theorem |5| solves the enumeration 



problem of counting non-intersecting lattice paths — by means of the determinant ( |2.3| ) 
totally. This is because \P {A ^ E)\ is easy to determine: Let A = (01,02) and E = (ei, 62) 



ILSE FISCHER 



be two points in Z^, such that A is located in the north-west of E (i.e., ei > ai and 02 > €2)- 
Then the number of lattice paths with steps in direction (1,0) and (0, —1) is 



\P{A^E) 



ei 



ai + 02 
ei - oi 



62 



ei 



ai + 02 

a2 — 62 



62 



(2.4) 



since each lattice paths corresponds to a choice of ei — ai steps in direction (1,0) out of the 
total number of ei — ai + 02 — 62 steps. 

Remark |l| shows that Theorem |5| can be useful in the enumeration of non-intersecting 
lattice paths with given initial and destination points. Since our enumeration problem of 
non-intersecting lattice paths not only involves fixed initial and destination points but also 
the additional condition about the fixed edge {x — l,y) — > (x, y). Theorem ^ seems useless in 
our situation at first glance. The following 'trick' remedies this matter: We add the following 
pair of initial and destination points: 



A 



a+l 



{x,y) and Ea+i = {x -l,y) (2.5) 

(see Figure ^). Then, as is not difficult to see, the number of rhombus tilings with a fixed 
rhombus 'at' {x, y) equals the number of families of non-intersecting lattice paths with initial 
points Ai and destination points Ei, i = 1, 2, . . . , a -|- 1. 

Now we are ready to apply Theorem^ to the points Ai and Ei, i = 1,2, ... ,a+l, defined in 
(p.l|), ( |2.2D and ( |2.5| ). In order to figure out what the determinant ( p. 3D actually gives, we have 
to find the permutations a for which there exist non-intersecting lattice paths connecting Ai 
to -Eo-(j) with l<z<a-|-l. It is quite easy to see that this is only accomplished by 
transpositions of the form (i, a + 1), where i ^ a -\- 1. The fact that all transpositions have 
the same sign — namely —1 — implies that the determinant in (|2.3| ) applied to our special 
points gives the number of families of non-intersecting lattice paths with initial points Ai and 
destination points Ei, 1 < i < a + 1, with negative sign. 

We use ( |2.4D to compute \P {Ai -^ Ej)\ in ( |2.3D , and finally obtain the following: 



Lemma 1. Let a,h and c he positive integers and {x, y) be an integer point such that < x < 
h + a — 1 and 1 <y < c + a — 1. Then the number of rhombus tilings of a hexagon with side 
lengths a,b,c,a,b,c which contain the fixed horizontal rhombus that corresponds to the point 
{x, y) in the bisection described above equals 



det 



1 < i < a 
( ^^^ ) 


J = 


= a + l 

-x+y\ 
-i+l) 


(c+x-y-i\ 
\ x-j J 








1 < i < a 



(2.6) 



In Lemma 1^ we refer to a correspondence between the integer points (x, y), < x < b+a—1 
and 1 < y < c + a — 1, and the horizontal rhombi we use for the rhombus tilings of a hexagon 
with side lengths a,b,c,a,b,c. This correspondence is implicitly given by the bijection between 
rhombus tilings and non-intersecting lattice paths described above. The following remark 
makes the bijection between points and rhombi more explicit. 

Remark 2. Again we consider a hexagon with side lengths a,b,c,a,b,c. We introduce the 
following oblique angled coordinate system: Its origin is located in one of the two vertices, 
where sides of lengths b and c meet, and the axes are induced by those two sides (see Figure ^). 
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The units are chosen such that the (Euchdean) side lengths of the considered hexagon are 
a,b,c,a,b,c in this coordinate system, too. (That is to say, the two triangles in Figure |^ 
with vertices in the origin form the unit 'square'). Thus, in this coordinate system, the 
points in Figure | have coordinates A[ = (1/2,9/2), A'^ = (3/2,11/2), A'^ = (5/2,13/2), 
E[ = (11/2,1/2), E'^ = (13/2,3/2) and E'^ = (15/2,5/2). Phrased differently, the coordinate 
system was chosen such that A[ = Ai + (1/2, 1/2) and E[ = Ei + (1/2, 1/2), where A[ and 
E[ are the initial and destination points of the paths in Figure |5|, and Ai and Ei are defined 
as in ( |2.1| ), ( |2.2D and (p.E|), 1 < i < a + 1. From this point of view, the family of paths in 
Figure |5| is the family of paths in Figure translated by (1/2, 1/2) and drawn in the oblique 
angled coordinate system. Furthermore, we see that the integer point (x, y) in the oblique 
angled coordinate system is just the lowest vertex of the corresponding rhombus under the 
aforementioned bijection between points and rhombi (see Figure ^ and Figure |^. 

3. From the determinant to a triple sum 

The aim of this section is the derivation of a triple sum that is equal to ( |2.6D and therefore 
gives the number of all rhombus tilings of a hexagon with side lengths a,b,c,a,b,c which contain 
a fixed rhombus with lowest vertex {x,y): 

Lemma 2. Let a, b and c be positive integers, and let {x, y) be an integer point such that 
< X < b + a—1 and 1 < y < c + a—1. Then the number of rhombus tilings of a hexagon with 
side lengths a,b,c,a,b,c which contain the fixed horizontal rhombus with lowest vertex {x,y) 
(in the oblique angled coordinate system; see Remark^ equals 

TT-rr -rr « + J + ^-l 1 (1) 



c + x — y + n — 2\fb — x + y + s — 1\ (^ + 1) 



n=l m=l s=l 



.1)«H 



1 J\ b-x + s-l J{b + c + l) 

^m - 1\ (c + l)n-l (6 + c + n)m- 



l) (n-1)! {m-n)\ J " ^^'^^ 

Outline of the proof of Lemma 12]. We derive the triple sum by starting from ( |2.6D . The 
matrix underlying the determinant in ( |2.6| ) has a 'homogeneous' definition, except for the last 
row and the last column. Our proof starts with some elementary row and column operations 
that transform the 'homogeneous' submatrix into a matrix of triangular form (see ( 3.10D ). 



Next we expand the determinant along the exceptional row and along the exceptional column 
(see ( 3.11| )). The result is a triple sum with a summand that involves another determinant. 



But this determinant is the determinant of the aforementioned triangular matrix with one 
row and one column deleted (which row and which column depends on the summation index 
of the triple sum, see (|3.13| )) and therefore we are able to compute it: The triangle property 
allows an expansion of this remaining determinant along the first (n — 1) columns and the 
last (a — m) rows, where m denotes the missing row and n denotes the missing column 
(see ( ^.14D ). Finally the remaining [m — n) x {m — n) determinants can be reduced to 
Vandermonde's determinant (see ( 3.15| )). 



Proof of Lemma |^ — The details. As mentioned before, we first describe some el- 
ementary row and column operations that transform the 'homogeneous' a x a submatrix 

--i+i) ) of ^^6 matrix underlying the determinant in ( |2.6D into an upper triangular 

/ l<i,j<a 

matrix. 
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We begin with some elementary column operations: By using the symmetry of the bino- 
mial coefficient (i.e., the identity (^") = („" J) in the last row of the matrix underlying the 
determinant in (p. 61) we observe that the determinant in (|2.6|) is equal to 



det 



1 < i < a 



( ^+' •) 

\C — t+J/ 



\c-y+j-lj 



j = a + l 

/b-x+y\ 
\y-i+lJ 



l<i<a 



(3.2) 



J i = a + l 



We add the (j — l)th column to the jth column, j = a, a — 1, ... ,2, in that order. The 
entries of the changed matrix read as 



b+c \ / b+c 
c-i + jj \c-i + j-l 



b + c+l 
c-i+j 



for i = 1,2, ... ,a and j = 2,3, . . . ,a, and 



c + X — y — l\ /c + x — y — l\ f c + x — y 

c-y + j-lj \c-y + j-2j \c-y + j-l 

for i = a + 1 and j = 2,3, ... ,a. The other entries do not change. Thus, the following 
determinant is equal to the determinant in (jsT 



det 

<i,jr<a+l 


i = i 
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\c-i+l) 


2 < i < a 

(b+c+l\ 
\c-i+j) 


j = 


= a + 1 

-i+iJ 




\ rr-V) 


/ c+x-y \ 
yc-y+j-l) 








1 < i < a 



(3.3) 



Next we repeat the procedure, i.e., we add the (j — l)th column to the jth. column, 
j = a, a — 1, ... ,3, in that order. Thus, we obtain that the following determinant is equal to 
the determinant in (t3.3|): 



det 

l<«j'<a+l 



i = i 
( ''^^ ) 

\c-i+lj 


J = 2 

(b+c+l\ 
\c-i+2j 


3<j<a 

(b+c+2\ 
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j = a + l 
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/c+x-y-l\ 
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(c+x-y\ 
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(c+x-y+l\ 
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1 <i < a 



We repeat this procedure of adding successive columns from right to left each time stop- 
ping one column earlier than before, as long as possible. This means that the procedure is 
performed a — 1 times including the two steps described in detail. Since the upper param- 
eter of the binomial entry increases by one every time it is involved and every entry in the 
jth column participates j — I times exactly, 1 < j < a, this procedure yields the following 
determinant: 
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det 



1 < i < « 

/b+c+j-l\ 
\ c—i+j ) 


j = a + l 

(b-x+y\ 
\y-i+lJ 


/c+x-y+]-2\ 
V c-y+j-1 1 






1 < j < a 



i = a 4- 1 



(3.4) 



Now we apply some elementary row operations to ( ^.41). In fact these row operations are 
analogous to the column operations we just applied to (^.S]). In order to do so, we first use 
the symmetry of the binomial coefficient for the first a rows of the determinant in ( |3.4D and 
observe that the determinant 



det 






1 < i < a 

(b+c+j-l\ 
\ b+i-1 J 



lc+x-y+]-2\ 



j = a + l 

I b~x+y \ 
\bM~x-l) 



l<i<a 



(3.5) 



J i=a+l 



is equal to the determinant in (BJJ). As announced, we now add the (i — l)th row of the 



determinant in (3.5) to the zth row, starting at i = a and stopping at i = 2. Next we do 
the same with the resulting determinant, starting at i = a but stopping at i = 3. After 
repeating this procedure o — 1 times we obtain that the following determinant is equal to the 
determinant in p? 



/ 



- det 



l<j<a 

b+i- 



(b+c+i+j-2\ 
-i-1 J 



/c+x-y+]-2\ 

y c-y+j-1 I 



j = a + l 

tb-x+y+i-l\ 
\ b+i-x-1 ) 



\ 



l<i<a 



(3.6) 



/ i = a 



Now we take the factor (-l)''+'=+*-i(6 + c + i- l)!/(6 + i - 1)! out of the ith row of the 
determinant in ( |3.6| ), z = l,2,...,a. This yields 

W 



TJ( ..b+c+i-i ib + c + i-l 
J-i^ ^ (6 + i-l)! 



X det 

l<«.i<a+l 



1 < i < a 

(c+j-1)! \b+c+i-l) 



\ 



tc+x-y+]-2\ 

^ c-y+j-l I 



j = a + l 



( -\\b+c+i-l (b+i-1)! /b-x+y+i-l\ 
V ^) (b+c+i-iy. \ b+i-x-1 J 



1 <i<a 

i = a + 1 

(3.7) 



Finally we want to apply the elementary row operations we just applied to the determinant 



in (3.5) once again. An analysis of the elementary row operations that we applied to the 
determinant in (p.5|) yields that the determinant in (p.6|) was produced from the determinant 
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in (3^) by replacing the ith row by 



T.('r])Ms) 



s=l 



s-1 



1 < i < a, where A{s) denotes the sth row of the determinant in (p]^). We perform this 
replacement of the entries in the determinant in (|3.7|) and obtain 



11(1 



b+c+i-l ' 



i=l 



(i-i)! ( i-j-i \ 

(c+j-iy. \b+c+i-l) 

— /,. I „ — 77-r7. T^ 



1 < j < a 



\^* (i-^\(_-\\b+c+s-lJHs^:2)^(b-x+y+s-l\ \l < i < a 
l^s=l\s-l!\ ^! (6+c+s-l)! I b+s-x-l ) \ 





j = a + l 



i = a + 1 

(3.8) 



The entry in the zth row and jih. column, 1 < i < a and 1 < j < a, develops from the 
corresponding entry in (|3.7D by using Vandermonde's summation formula: 



E 



i\ (i-1)! 



(i-1)! 



l) {c + j - l)\\b + c + s - l) {c + 3 - l)\\h + c + i - l) ' 



i-1 



Next we apply the elementary identity 



— n + A; — 1 
k 



(3.9) 



to this entry in the ith row and jth column of the matrix underlying the determinant in (|37 
1 < i < a and 1 < j < a, and then use the symmetry of the binomial coefficient. Finally we 
take the factor (— l)''+'^+*^i out of the ith row, 1 < i < a, and obtain that the expression in 
(^^) is equal to 



n 



(6 + c + i-l)! 



L {b + i 



l)\ 



X det 



1 < j < a 

0-1)! /fc+c+j-n 

(c+j-l)! I j-i ) 



tc+x-y+]-2\ 



El ^*^i 



j = a + l 
-^Y^s jb+s-iy. ^b-x+y+s-l\ \l<i<a 



(b+c+s~iy. \ b+s-x-l 







i = a + 1 

(3.10) 



With pleasure we discover that the a x a submatrix induced by the first a rows and first 
a columns of the matrix underlying the determinant in ( p. 10 ) is an upper triangular matrix. 
This is because the binomial coefficient (^) is defined to be zero if /c < for any indeterminante 
a. Therefore the determinant is of that form we were looking for at the beginning of the 
proof. 
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Next we expand ( |3.10| ) along the last row and then along the last column, to obtain 

n (b + c + i-l)! A ^.a+i+n 
(b + i-l)\ ^^ ' \ c-y + n-1 

i=l ^ ' n=l ^ 



{h + c + i -JO! y^/_;LNa+i+n fc + x-y + n-2 

n=l ^ 

ir'"E(7:i')(-i) 



x^/ .^a-^^.^ fm-l\ _^ (6 + s-l)! (h-x + y + s-l 



^^' ' -^^Vs-ly (6 + c + s-l)!V 6 + s-x-l 

ni=l s=l 

i<ij<a V(c+j-i)!v j-« yy 

Then we take the factor (j — l)!/(c + j — 1)! out of the jth column of the remaining 
determinant. This, together with some other manipulations, gives 



n 

i=l 



+ i-l)\{c + i-l)\ 

a a m 

EEE(-ir 



„, c + x — v + n — 2\ fb — x + y + s — 1\ fm — 1 

n=l m=l s=l 



(» + "-')! ('= + "-^)! det (■f' + '^/^-'n. (3A2) 



(6 + c + s-l)! (n-1)! i<jj<a \\ j-i 



Now we have to compute 



det (('^-'^^ {-')]. (3.13) 

i<ij<a VV j-^ yy 

This is a determinant of an upper triangular matrix, where the mth row and nth column 
was deleted. Therefore it is only different from zero \i n < m. So, let us assume n < m. 
Expansion of the determinant along the first (n — 1) columns and along the last (a — m) rows 
yields 

det (C' + =, + ^:-l))= det fC + ' + f)) (3.14) 

n+l<j<m\\ J—t J J n<j<m-l \\J — Z + ly y 

n<i<.m—l n<.i<m—l 

for the determinant in ( |3.13| ). 

We are going to reduce this determinant to Vandermonde's determinant. In order to do 
so, we take (6 + c + j)!/(j — n + \)\ out of the jth column, 1 < j < a, and 1/(6 + c + i — 1)! 
out of the ith row, 1 < i < a. Thus, we obtain that the determinant in ( |3.14|) is equal to 



m—n 



{b + c + n)m-nT[-7 det ((j - i + 2)i„i) . (3.15) 

J- l<]<m—n 
■'~ l<i<m—n 

The entries of this determinant are monic polynomials in j of degree i — 1, where i and j 
denote as usual the index of the row and the column of the entry. It is now straightforward 
to reduce this determinant by appropriate row operations to Vandermonde's determinant, 

(3.16) 









m—n—1 


^^det ((j- 


-i + 2),-i) = 


= det (/-I) = 

l<j<m-n ^ ^ 


= ^ i\ 


l<j<m—n 




l<j<m—n 


i—l 
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The last equation was obtained by applying the well-known formula for Vandernionde's de- 
terminant to deti<j<m-n,i<i<m--n (f^^)- 

Using ( |3.16| ) in ( p. 15 ), we obtain that the determinant in ( p.l3|) equals 



m—n , m—n—1 



{b+c+n)m-n n ~ n 



+ c + n) 



j=l •> i=\ 



(m — n)! 



Accordingly, we replace the determinant in ( 3.12 ) by its value (6 + c + n)m-n/(^T2 — n)!. Thus, 
we obtain the following triple sum for the determinant in ( |2.6| ) : 



n 



(b + c + i-l)!(i-l)! 
+ i-l)!(c + i-l)! 



a a m 
n=l m,=l s=l 



(-1)"' 



c + x — y + n — 2\/6 — x + y + s — l\ (m — 1 



X — \ j \ 6 — x + s — 1 j \ s — \ 

(6 + 1)^-1 (c+ l)n-i (6 + C + n)m-n 
3 + c+l)s-i (n — 1)! {ra — n)\ 



Therefore Lemma ^ is finally proved since 



n 

i=2 



' + c + i-l)!(z-l)! 



a h c 

nnn 



i + j + /c - 1 \ (1) 



(6 + i-l)!(c + i-l)! lllllllz + j + /c-2 I (6+1 



4. Evaluation of the determinants 

In this section we compute the determinant from Lemma || in case that the fixed rhombus 
is placed in the centre (see Lemma p| and in case that the lowest vertex of the fixed rhombus 
is placed in the centre (see Lemma^) . The combination of these two lemmas and Lemma IT] 
then establish Theorem IT] and Theorem 0. 

In order to do so, we first have to figure out which integer point (x, y) corresponds to 
the central rhombus, respectively the 'almost central' rhombus above the centre, via the 
bijection described in Section Q. In the oblique angled coordinate system introduced in 
Remark Q the point ((a + 6)/2, (a + c)/2) is the centre of the considered hexagon. Because 
of that, ((a + 6)/2, [a ^ c— l)/2) is the lowest vertex of the central rhombus and therefore 
the integer point that corresponds to the central rhombus (see Remark ^). Accordingly, 
((a + 6)/2, (a + c)/2) is the lowest vertex of the 'almost central' rhombus above the centre. 

The evaluations of both determinants, namely the determinant corresponding to the case 
that the fixed rhombus is placed in the centre, and the determinant corresponding to the 
case that the fixed rhombus is placed next to the centre on the other hand, are quite similar. 
Therefore we concentrate on the first case, see the following lemma. Lemma Q at the end of 
this section is devoted to the second case, but there I only explain the major differences to 
the first case. 



Lemma 3. Let a,b,c be integers. Then the determinant 
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det 

l<«,i<a+l 



/ 


1 < i < a 
f ^ + ^ ^ 

\c-i+ j) 


j = a + l 
( '-^ \ 

Vf-i + ^y 


\ 


/ b+c-l \ 






1 < j < a 



i = a + 1 



is equal to 



(a b c 
nnn 
i=ij=ik=i 

(a-l)/2 

X E 

A:=0 



i+j + k-l\ (1), 



1 



b+c-l \ /a+ b+c-2 

2 \oa-l 

b-1 l\ 6-1 '^ 



i+j + /c-2 / (fe)e(6 + C+l)a_iV V y V ^ 



c+l\ /l + 6 + c\ /c + 2fc + 2 



X 

in case that a is odd, and 
, 6 



6 + c + 2/i: + 3 



(a-2fc-l)/2 

(2)(a-2fc-l)/2 
(a-2fe-l)/2 (l)(a-2fc-l)/2_ 



i + i + A;-l \ (1), 



^JlUi'+^ + ^-^l ib)c{b + C+l)a-l 
(a-2)/2 

E 



c + 2\ /l + b + c\ /c + 2k + 3 



k=0 



in case that a is even. 



b+c-l \ / a+b+c-l 
2 
b 
2 ^ \ 2 



(a-2fc-2)/2 



2 |2"-2 



6 + c + 2A; + 3 



(2)(a-2A;-2)/2 
(a-2fe-2)/2 (l)(a-2A;-2)/2 



Outline of the proof of Lemma |3|. The following procedure is used for computing the 
determinant: The determinant depends on the side lengths o, b and c of the hexagon. More 
precisely: The dimension of the matrix underlying the determinant is a + 1, and the entries 
are certain binomial coefficients that depend on b, c and a. First, in Step 1, we reduce the 
problem to the computation of a determinant with entries that are polynomials in b and c 
if we fix a, see (|4.6|), and therefore a determinant that is itself a polynomial in those two 
variables. The comparison of the determinant in Lemma ^ and the determinant in ( |4.6| ), the 
latter will be denoted by det{Da(b,c)), yields that we have to show that the determinant 
det{Da{b,c)) is a product of certain linear factors and an irreducible polynomial (over Z) in 
b and c, see ( [4.7D , respectively ( |4.8| ). 

In Step 2 of our proof we show that every linear factor on the right-hand side of (4.7), 
respectively (4^), is indeed a linear factor of det(Da(^) c)). I explain the used procedure by 
an example: We have to show, e.g., that the linear factor (c + 1) is a factor of det{Da{b, c)) 
in case that a is even (see ([4.8|)). It is a fundamental algebraic fact that, in order to show 
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that, it suffices to sfiow tliat det(Z?a(6, — 1)) = if a is even. (Here we use tlie fact that 
det{Da{b,c)) is a polynomial in c if we fix a). Clearly, a determinant of a matrix is equal 
to zero if and only if there exists a linear combination of rows, or, equivalently, of columns. 
Therefore we find a linear combination of rows of Da{b, —1) which vanishes, and then prove 
it, see ( [4.13| ). Proving in this case means to establish hyper geometric identities Q. 

At this point it is worth mentioning that the procedure described so far gives a complete 
solution for determinants that factorise completely into linear factors. 

In Step 3 of the proof we finally compute the irreducible polynomial, which will be denoted 
by Pa{b,c). In order to do so, we first look for special values of b, where Pa(b,c) is 'nice'. 
Indeed, we discovered that Pa(b,c) factors completely into linear factors for b = —c — k, 
/c = 1, 3, . . . , 2[(a - 1)/2J + 1. We work out these evaluations of Paib, c) in ( pOD and (|^), 
and subsequently prove them by making use of the triple sum derived in Section 0. As it 
turns out, the degree of Pa{b, c) as a polynomial in b is exactly [(a — l)/2j . Thus, the above 
[(a — l)/2j + 1 evaluations suffice to compute Pa{b,c) by using Lagrange interpolation, see 
(PI and iH). 

Proof of Lemma ^ — The details. 

Step 1: From our determinant to a determinant with polynomial entries. 

In Section g it was finally shown that the number of rhombus tilings which contain a 
fixed horizontal rhombus with lowest vertex (x, y) is given by the following (a + 1) x (a + 1) 
determinant: 

1 < j < a j = a + 1 



det 



\C — l+]J 



tc+x~y—V 
\ x-j 



) 



/b-x+y\ 
\y-i+lJ 



l<i<a 



(4.1) 



With displeasure we observe that the matrix underlying this determinant has an excep- 
tional row and an exceptional column, namely the (a + l)th in both cases. In the following 
I describe some elementary row and column operations which lead to a matrix having the 
same determinant but a 'homogeneous' definition: 

1. First we want the entries in the (a-l-l)th row to be zero — except for the first entry. Since 
{'"'^Ij''^) is the jth entry in the (a + l)th row, we have to subtract {"'^IZ^'^) / ("t-'Xi^) = 
{x — j + l)/(c — y + j — 1) times the (j — l)th column from the jth, starting at j = a and 
stopping at j = 2. We obtain a new matrix with the desired behaviour in the (a + l)th row. 
The {i,j) entry of this new matrix is 



b + c 
c-i+j 



X 



J + 1 



y + j-l\c 



b + c 
i+j- 



1 



(c+x-y-l\ 



a + 1 and j = 1, and it is for i = a + 1 and j ^ 1. 



t,j = 2,3,...,a,itisn-_r')foT 
The other entries do not change. 

2. Next we want to do the same for the (a-|- l)th column: Accordingly, since ( ~J^Tf) is the 

ith entry in the (a + l)th column, we subtract (^"^+j')/(yl^+|) = (y - i + 2)/{b + i-x-l) 
times the (j — l)th row from the jth row, again starting at the bottom — that is to say, 



^At this point it is worth mentioning that I used C. Krattenthaler's Mathematica package HYP 
most of the hypergeometric identities within this paper. 



to handle 
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j = a, a — 1, ... ,2. Therefore the (i, j) entry of our new matrix is given as fohows: 

b + c \ X — j + 1 / b + c 

.c-i+jj c-y + j-l\c-i + j-l, 
y — i + 2 ( ( 6 + c ^ x — j + 1 / 6 + c 



6 + i 



1 



1)+J 



y + j - 1 Vc- (i- 1) + j - 1 



(4.2) 



if i, J = 2, 3 . . . , a, it is ( ~^ ^) for j = a + 1 and i = 1, and it is for j = a + 1 and i ^ \. 
Again, the other entries do not change. 

3. In summary, we obtain the following determinant. 



/ 



det 



V( 



* 

x-1 ) 







/b-x+y 









)\ 



/ 



where * describes an entry which is of no interest for us, and where aij denotes the expression 
in (|4.2| ). This form suggests an expansion of the determinant along the last row and along 
the last column. Therefore (|4.l| ) is equal to 



c + X — y 
X — 1 



l\ (h — X + y 

y 



det (aij). 

2<t,j<a 



(4.3) 



We have reduced our problem to computing det2<ij<a{aij) and therefore the job to com- 
pute the determinant of a 'homogeneous' matrix, even if the entries are now more complex. 
The advantage is that we can now take several factors out of the determinant so that the 
remaining entries are polynomials: A simple calculation shows that 

(6 + c)!(6 + i-j + 2)j-2(c -i + j + 2)„_ 
"""^ {b + i-iy.{c-j + a + l)!(c -y + j-l){b + i 

where H is the polynomial 



X 



1) 



H = {b + i-j + l)(c - i + j + l){c - y + j - l){b + i - X 
— {c — i + j){c — i + j + l){x — j + l){b + i — X 
-{b + i-j){b + i-j + l)iy-i + 2){c-y+j 



X H, 



1) 
1) 
1) 



(4.4) 



+{b + i-j + l)(c -i + j + l)(x -j + l)iy-i + 2). 



We take {b+c)\/ {{b+i-l)\{b+i- x-l)) out of the zth row, and l/((c- j + a-l)!(c-2/+j-l)) 
out of the jth column of the matrix {aij)2<ij<a- This gives 



det (aij) 

2<t,j<a 



n 

\i=2 



+ c)\ 



' + i- l)!(c -i + a + l)!(c - y + i - l)(b + i - x - 1) I 
X det {{b + i-j + 2)j.2{c-i + j + 2)a-j-H). 

2<i,j<a 



(4.5) 



Again we have reduced our problem to the evaluation of another 'homogeneous' determi- 
nant, namely 



det {{b + i-j + 2)j_2(c -i + j + 2)a-j ■ H) 

2<i,j<a 



(4.6) 
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But this one has the pleasing property to be a polynomial in b and c if we fix a (since the 
entries of the underlying matrix are). In the following, Da{b, c) denotes the matrix underlying 
the determinant in ( [4.6| ) evaluated at x = (a + 6)/2 and y = {a + c — l)/2. 

When I computed this determinant for small values of a, I was led to the following conjec- 
ture 0: There holds 

.<t^. P.(^ c)) = (ll (1 + 5 + c),_,^ (^f[i^ - 1)!^ 2-^ 
b + lV (c+2\ fl + b + c 

/ {a-l)/2 

k 



2 

(a-l)/2 

E 



A:=0 



(a-l)/2 \ ^ / (a-l)/2 

l + b + c\ fc + 2k + 2\ 

k\ 2 / (a_2fe_l)/2 

6 + C+2A; + 3\ (2)(a-2fc-l)/2 



X 



{a-2fe-l)/2 (l)(a-2A:-l)/2 



(4.7) 



if a is odd, and 



det {Da{b,c)) 

2<i,j<a 



lia + b + c),} (fl{i-l)!\2^ 
/c+1 



U=2 



6 + 2 



l + b + c 



(a-2)/2 

X E 

fc=0 



(a-2)/2 \ ^ / a/2 V ^ /a/2 

c + 2\ /l + 6 + c\ /c + 2A; + 3 



6 + c + 2A; + 3 



(a-2fe-2)/2 

i2){a-2k-2)/2 
{a-2k-2)/2 (l)(a-2A:-2)/2 



(4.8) 



if a is even. If we remember the assertion in Lemma y and the factors we have taken out of the 
determinant in ( ^ ) and (^^), and if we then specialise x = {a + b)/2 and y = (a + c— l)/2, a 
tedious but straightforward calculation yields that in order to complete a proof of Lemma |3| 
it remains to show ( [4.7| ) and ( [4.8]) . 

As described in the outline of the proof, we first prove that the claimed linear factors, 
which are 



n(i+^+c)._iV 



b+l\ fc + 2 

(a-l)/2 



(a-l)/2 



l + b + c 



•,2a-2 



(4.9) 



(a-l)/2 



^In fact, only the linear factors in front of the sums in (4.7) and (4.8) can immediately be guessed, by 
computing det{Da{b,c)) for a = 2,3, ...,8. In order to work out an conjecture concerning the two sums 
themselves I first computed them for a = 2, 3, ... , 23. On the basis of these data, I worked out (guesses for) 
'nice' evaluations of these polynomials at special values of 6, for all values of a. (See the Outline of the 
proof of Lemma H and Step 3 of the proof of Lemma M for more details, and, in particular, for an explicit 
hsting of these special values of h) . There by I was helped by Krattenthaler's Mathematica 'guessing machine ' 
RATE (which is available via Internet at ittp: //radon. mat .univie . ac. at /Peopl e/k ratt / rat e/rate. html; 
see also ISl, Appendix A]). Once I had got thus far, I computed the polynomials in (4.7) and ( }4.8[ ) by Lagrange 
interpolation. 
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if a is odd, and, respectively, 

(llii+b+c),_Ab 



6 + 2 



c+1 



1 + b + c 



-,2a~2 



(4.10) 



^j=2 / \ - / {a-2)/2 \ ^ / a/2 V ^ /a/2 

if a is even, divide det(Z?a(^, c)) (see Step 2). In Step 3 we finally calculate the remaining 
irreducible polynomial in b and c (over Z), which reads 




n( 

\i=2 



if a is odd, and respectively, 

\i=2 

if a is even. 



c+1 



1 + 6 + c 



6 + c + 2A; + 3 



c + 2A; + 2 



(a-2fc-l 



(a-2fc-l)/2 

(2)(a-2fc-l)/2 
)/2 (l)(a-2fc-l)/2 



(4.11) 




c + 2 



1 + 6 + c 



+ c + 2A; + 3 



c + 2A; + 3 



(a-2fc- 



(a-2fe-2)/2 

(2)(a-2fe-2)/2 
-2)/2 (l)(a-2fc-2)/2 



(4.12) 



Step 2: The linear factors in (4.9), respectively in ( ^.1C1| ), divide det{Da{b,c)) as 
a polynomial in b and c. 

Essentially we have four different types of linear factors: 

1. Factors of the form (c + k): If a is odd, k = 2,4, ... ,o— 1. Otherwise A: = l,3,...,a — 1. 

2. Factors of the form (b+k)'^: If a is odd. A; = 1, 3, . . . , a — 2. Otherwise A; = 2, 4, . . . , a — 2. 
The factor b occurs once if a even. 

3. Factors of the form (6 + c + k)''-'' for A; = 1, 3, . . . , 2[(a - l)/2] - 1. 

4. Factors of the form (6 + c + kf"''-'^ for A; = 2, 4, . . . , 2[(a - l)/2j . 

re 1. — The factors of the form (c + A;) divide det{Da{b, c)): In order to show that (c + A;) 
is a linear factor of det(-Da(^5 c)), where 1 < A; < a and k ^ a (mod 2), it suffices to show that 
det(-Da(6, —k)) = for those special A;'s. Therefore we search for linearly dependent rows or 
columns in the matrices Da{b, —k). 

By computer experiments PI , I found the following non-trivial linear combination of rows: 

a-fe+l (—\\ 



E 

i=(a-A;+3)/2 

for J = 2, 3, 



+ l)a-k+l- 



'a+k~l+2i \ 

2 I a-k+l-i 



(l)a- 



■fc+1- 



'fe-a+2^-2^ 



-Da{b,-k)(^ij) = 0, 



(4.13) 



i—k+l—i 



,a. 



In order to complete the proof that (c + A;) divides the determinant det{Da{b, c)), we just 
have to prove this identity. Gosper's algorithm for hypergeometric sums (which are sums, 
where the quotient of two successive summands is a rational function in the summation 
index) helps us to recognize that this is actually a telescoping sum: First we examine the 



''Again, I computed these linear combinations for small values of a, i.e., I solved the follow- 
ing system of linear equations in Ci(a,b,k): ^"^2 Ci(a,6, fc)Da(b, — fc)(ij) — 0, j — 2,3, ... ,a, for 
a = 2, 3, . . . , 15. I was led to an conjecture for the coefficient Ci{a,b,k) for all values of a 
by using Krattenthaler's Mathematica 'guessing machine' RAT E (which is available via Internet at 
tittp: //radon. mat .univie .ac . at /People/kratt /rate/rate. html). 
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factor ((—a + k — 1 + 2i)/2)a-k+i-i in the linear combination. I claim that this factor is 
equal to zero if i < (a — A; + 3)/2 and because of that we can omit the lower bound on the 
summation index on the left-hand side of ( [4.13|) : In this case (—a + A; — 1 + 2i)/2 < and 
(a — k — l)/2 > 0. Because of that, and since (—a + A; — 1 + 2i)/2 is an integer, one of the 
factors of the Pochhammer symbol 

-a + A;-l + 2i\ _ / -a + k - I + 2i\ / -a + k - 1 + 2i \ fa-k-l 

2 )a-k+l-^ ^ V 2 ) \ 2 + V • ■ ■ V 2 

is equal to zero. If we then reverse the order of summation, we obtain that we must show 
the following identity: 

^ (-iyQ-f + |)J-a-fe + fc).+,^2 ^, 

/ . TTT T, a h , i\ -" \x={a+b)/2,y=(a+c-l)/2,c=-k,i-,a-k+l-i - U. 

(4.14) 



Let f{i) denote the summand of the sum in (4.14). The magnificent algorithm due to Gosper 
decides whether there exists a hypergeometric g{i) with 

g{i + l)-g{i)=f{i). 

And in case of its existence the algorithm also computes g{i). If such an g[i) was found we 
would have 



oo 

E 

i=0 i=0 



f{i) = E ^9{i + 1) - 9{i)) = lim {g{i) - ^(0)) . 



(limj_»oo5'(i) exists in our cases, since f{i) = for all but a finite number of i and therefore 
g{i) is finally constant.) 

A computer implementation of Gosper's algorithm p^ prints out 



(-l)n-fl - h + fc),+j_2 (^ - § + f ), (-1 - b + fc)(i - j + fc) 
2(l)_i_,+i+, (l-f-| + A;)^_^ 

to be a suitable g{i) for our /(i). One may check the identity g{i + l)—g{i) = f{i) by dividing 
the left-hand side by the right-hand side and simplifying the resulting rational function to 1. 



This implies ^ f{i) = 0, i.e., the truth of ( 4.14 ), since g{Qi) = 0, caused by the factor 

j=0 

l/(l)_i_a+j, and g{i) = for i > (1 -|- a — k)/2, caused by ((1 -|- a — A; — 2i)/2)j (again we 
use the fact that 1 < A: < a and k ^ a (mod 2)). Thus we have proved that c -|- A: is a factor 
ofdet(i^a(^c)). 

re 2. — The factors of the form (b + k) divide det{Da{b,c)): Next we have to prove that 
(6 -|- A;)^ is a factor of det(Da(6, c)) if < A; < a and k = a (mod 2), and, in addition, that b 
is a factor of det{Da{b, c)) if a is even. 

In case that A; = a — 2, this is easy: If we examine the matrix Da{—a + 2, c) carefully, we 
observe that the (a — l)th and ath row vanish. (For j > 2 this is caused by {b + i — j + 2)j_2, 
and for j = 2 by the nasty polynomial H.) Therefore (6 -|- a — 2) is a factor of the (a — l)th 
row and a factor of the ath row of Da {b,c). Thus, we have proved that (6 -|- a — 2) ^ is a factor 
of the determinant det{Da{b,c)). 

In case that k < a — 2 we use the same method as for the factors of form (c -|- k) above. 
But since we are now dealing with factors of higher multiplicity, it is not enough to find just 
one linear combination of rows or columns of Da{—k, c). In fact we have to find two linearly 
independent combinations if we want to prove that (6 -|- k)'^ is a factor. 
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We claim that 



i=k+2 UA-fc-1 I 2 )i-k~2 



)i-k- 



xpi_k-i{c + a-k-i + l)Da{-k,c)(^ij) = (4.15) 



ifO<A;<a — 2 and k = a (mod 2) , and that 



(.+^2)/2 (_iy-/c(, + a-i + 2),_,_i (a±Al_2l±4; 
i=A;+3 



i-k-l 



/i\ ( c+a-2i+3 \^ 

yi^Ji-k-l \ 2 Ji-k-2 



X pj_fc_2(c + a - A; - i)Da{-k, c)(jj) = -Da{-k, c)(^k+i,j) (4.16) 



ifO</i;<a — 2 and k = a (mod 2) are such hnearly independent hnear combinations of 
rows, Pn{c) being the sequence of polynomials given by 



ra-l 



Pn{c) = J2 



h=0 



1 + c — n 



n—h—l 



1 + c-2h + n 



(4.17) 



Notice that the exceptional factor b of det{Da{b, c)) for a is even, is included in the first linear 
combination ( 4.15 ). 

In order to see that the two linear combinations in ( [4.15| ) and ([4.16| ) are linearly independent 
we remark that the first linear combination ( [4.15| ) involves the {k + l)th row of Da{—k,c) 
whereas the second linear combination ( [4.16| ) does not (here we use k < a — 2). 



It remains to show the two hypergeometric identities (4.15) and ( 4.16| ). The situation is a 
bit more complicated this time compared to the hypergeometric identity (|4.13| ) that proved 
that (c + A;) divides det (Z?a (6, c)): Actually these two identities are double sum identities, 
since they involve the polynomials Pn{c). Luckily the proofs of the two identities are quite 
similar. We start with ( 4.15| ). 
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To begin with, we split the sum into four sniaher sums — each one corresponding to one 



of the four summands of H (according to the representation of H in (4.4)): 



<'"+»+'>/= (-l)*-'-2(c + a -i + 2), ■ '■'+'-M+1' 



E 

i=k+2 



li-k-2 



i-k-2 



\.^)t-k-l \ 2 Ji~k-2 



E 



X Pi-k-i{c + a-k-i + l)Da{-k, c)(jj) = 

Pi^k-iic + a - k - i + 1) 



(a+^2)/2 (_i).-.-2(^ + „ _ . + 2),_,_2 {B±h^l±±)^_^_^ 



(1). 



'c+a-2J+3^2 



i=k+2 \^)i-k-l V 2 n-k-2 

x{-k + i-j + l)j-_i(c - i + j + l)a-j+i((-l - a + c + 2i)/2)((-2 - a + 2i - A;)/2) 



(„+,+2)/2 (_i).-.-2(^ + „ _ . + 2),_,_2 (^^±^^2^),_,_2 



+ E 



/c+a-2i+3\2 
fc+2 V-LjJ-fc-l I, 2 li-k-2 



(l)i-fc- 



Pj_fc_i(c + a- A;- i + 1) 



-A: + i - j + 2)j_2(c - i + j)a-j+2{i-2 -a + 2i- k)/2){{2 + a - 2j - k)/2) 



+ E — 



( a+k-2i+4: \ 
h-k-2 \ 2 )i-k-2 



(1). 



c+a-2i+'i 



•A\2 



Pi-k-i{c + a- k- i + 1) 



=fc+2 K^)i-k-l \ 2 Ji-k-2 

x{-k + i- j)j{c -i + j + 2)a-j{{3 + a + c- 2i)/2)((-l -a + c + 2j)/2) 



^-+^+^y\-iy-k-2{c + a~i + 2), ■ ^^-tfc^2^+4^ 



+ E 



)i-k-2 



i-k-2 



[^')i-k-\ 



'c+a-2J+3\2 



Pi-k-i{c + a- k-i + l) 



i=k+2 \^Ji-k-l \ 2 Ji-k-2 

x{-k + i-j + l)j-_i(c -i + j + l)a-j+i((3 + a + c- 2i)/2)i{2 + a - 2j - k)/2). (4.18) 



Next we examine the sum of the first and the third sum of the right-hand side in (4.1^ 
After some simphfications and interchange of summations one gets 



{a-A;-2)/2 

E 

h=0 



{a+k+2)/2 (_l _ a , c 

^y-j+l { 2 2 



E 



+ -2+j){k--2--2+k). 



+ k), 



3 a c 

i=h+k+2 V2 2 2 ' 'Vh 

(1 - f + l)-l+i-k (1 - a - C + k)-l+a+^-j-k 



Wi-j-k (2 2 2 + ^j- 



'2+i-A: 



(a-fc-2)/2 

E 

h=0 



\a+k+2)/2 
i=h+k+2 



' '^ + i+j)(l-f-i + fc). 



2 2 



.2 2 2 "•" '^M 



1 - § + I) -2+^-fc (1 - Q - C + fc)_2+a+i-j- 
(l)-l+i-i-fc (2 ~ f ~ f + ^)-3+i-fc 



Notice that these two double sums are quite similar. If we combine the exterior sums, take 
some common factors out of the two inner sum and shift the index of the first inner sum by 
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one, we obtain 

{a-k-2)/2 



\2 2 2 ^ '^Jh 



h=0 



\a+k+A)/2 /-, a j_ k\ (-i „ ^ _i_ i-^ 

1 - 2 + 2 J -2+i-k (1 - a - C + fcj„2+a+i-i-fc 



E 

i=h+k+3 



l+i-j-k\2 2 2+^)-3+J- 



(1) 
(a+fc+2)/2 /, 

E - 

j=/i+fc+2 






2/-2+j-fc 



(1 - a - C + k)^2+a+i-j-k 



(1) 



-1+i-j-k [2 2 2 ~^ ^) ^ 



-3+i-k 



This makes us discover, that the two inner sums nearly cancel out each other, because the 
two summands are exactly the same. Therefore the latter expression simplifies further to 



(a-A:-2)/2 
h=0 



i — Ql J- £ J- n\ (X — °L — £. 
2 2 ^ 2 ^ -J I \2 2 2 



+ k). 



3 a c 

2 2 2 



+ k), 



(l-§ + l)h(l-°-C+^)a+fc- 



)l+h-j (f-I -§ + ^)_i_ 



(l)l+/i-j \2 2 2 ' '"J-l+h 

Thus we have already reduced the sum of the two double sums to the following single sum 

h (i-f-§+fc)Jl)m-. 

(4.19) 

Analogously one can show that the sum of the second und the fourth sum in ( [4.18| ) , which 
reads 



(a-A:-2)/2 (a+fc+2)/2 

E E 

h=Q i=2+h+k 



k^ ^1 — 9l — £ -i- lA 
2 2 2 ^ '^)h 



■K-z]/z ^a+K+zj/z (^ \ a _ ■ _ k\ I 

V^ V^ l..-l+a-jV_^J2^_!__2n 

l^ Z^ ^ ' ('3_a_c , r.\ 

\2 2 2 ^ '^Jh 



X 



:i-f+i)-i+.-fe(i-«-c+A;u,_,_, 



(l)l+i-i-fc (2 " '^ 



2 2 + ^) -2+i-k 
-K-zj/z (a+K+z)/z (^ \ a _ ■ _ k\ ( 



(a-A:-2)/2 (a+fc+2)/2 

E E 

h=Q i=2+h+k 



k\ /I a c 

2 2 2 



+ k), 



X 



(^ f + 2)-2+J-A: (-*- ^ '^ + k)-l+a+i-j-k 



(1) 



i-j'k (2 2 2 +^)-3+i-fc 



2 2 ' 'V-3+i-fc 

after interchange of summation and some cancellations, simplifies to 

h (iW/.-.(i-f-§+^). 

(4.20) 

Now it remains to show that the two (single) sums ( 4.1S| ) and ( 4.2[1| ) sum up to zero. 
Because of the factor l/(l)i+/j_j, we are able to change the lower bound on the summation 
index of the sum in ( 4.1S| ) to j — 1, since j > 2 and l/(l)n = if n is a negative integer. 
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Analogously we change the lower bound of the summation index of the sum in ( [4.20 ) to j — 2 
according to the factor l/{l)2+h-j- Therefore the sum of ( 4.19| ) and ( [4. 201 ) is 



lac 



2 2 2 



+ k 



(a-fc-2)/2 

E 

h=j-l 



(_^)a^,lzi^i_+§+j')(l-Ul).(l-«-C + fc)a+h-, 



.2 - I - i + ^)/i(l)i+/i-i 



"''J2'\_,^a-,{2±i^l_- I) (1 - f + I). (1 - a -c+k),^^ 

h=j-2 



+h~j 



)2+h-j (2 ~ f ~ § + ^)ft 



Using the standard hyper geometric notation 



r-' s 



ai,. . . ,ar 
bi,.. 



\z 



E 

fc=0 



(ai)fc . . . {ar)k k 
k\{bi)k...{bs)k 



and, after canceling the factor (^^ — ^ — ^ -\- k^, it remains to show that 



(1 



2 "•" 2' 



-1+i 



(1 



c + k)^ 



-l+a 



+ k)-i+j 



2F1 



" +j + '' 



1' a 



2' 



-C + k 



-2-f'l 



l+j + k 



1 + 1+J) 



2 ' J ' 2' 

a c 

2 2 2 



-C + k 



^ " " +j+k 



+ j + /c) =0. (4.21) 



The factor ((2 — a + A;)/2)_i+j tells us that this is true at least for j > (—2 + a — k)/2. 
Otherwise, we use Vandermonde's summation formula (see [|^, (1.7.7), Appendix (III. 4)]) 



2F1 



a, —n 



;i 



(c - a)r 



(4.22) 



which is valid if n is a nonnegative integer. Namely, if we apply this formula to the left-hand 
side of ( 4.21D with n = {a — 2j — k)/2, respectively n = {2 + a — 2j — k)/2 (here we use 
j < (—2 + a — k)/2), we see that the two 2-^i^series in ( 4.21]) sum up to zero. Thus, we have 
proved the first linear combination, ( 4.15| ), for the factors of form [b + k). 



The proof of the second linear combination ( 4.16| ) is quite similar: Again we split the double 
sum into four smaller sums according to the summands of the polynomial H. After some 
simplifications and interchange of summation we have reduced our problem to the following 
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hypergeometric identity: 

(-4+a-fc)/2 (2+a+fc)/2 



E 

i=3+h+k 



E 

h=0 

(-4+a-fc)/2 (2+a+A;)/2 

E 

h=0 i=3+h+k 



(-1) 



1+a+h-j 



+ i + j)i (^ + |)i-fc {-a-c + k)a+i-j-k 



+ 



+ E 

h=0 
(-4+a-fc)/2 (2+a+fc)/2 

E E 

h=0 i=3+h+k 



^)i—j—k 
1) 



3 I g I c 

'2 "•" 2 "•" 2 



h - k)i+h ( 



+ h + k)_2-h+i-k 



l+a+h—j /I _|_ a 



■ J - |)l (^ + l)i~k {-a-C + k)i+a+i-j' 



(l)l+j-j-fc 



'2 "•" 2 "•" 2 



^)l+/i (211 



+ h + k)_2-h 



+i-k 



f + § + j)l (^ + 2)-l+«-fc ("'^ — C + k)^i+a+i-j~k 



(1) 



-1+i-j-k 



3 I a I c 



/l - A;)i+h ( 



, 3 

l+h{2 



+ h + k)^3-h+i~k 



+ 



(-4+a-fc)/2 (2+a+fc)/2 /_,Nfi 



E 



(i + f 



i - |)i (^ + |)-i+i-fc (-a - c + k) 



a+i—j—k 






h=0 



{l)i-j-k (- 



i=3+h+k V-^^-J-f^ V '2 "^ 2 "^ 2 



/i - k)i+h (|-f-§ + /i + A:)_3_/j+i_fc 

= -i?„(-fc,c)(,,+i,,). (4.23) 
As before, the sum of the first and the third summand on the left-hand side of ( [4. 231) 



simphfies to a single sum: 

,a+fe-, (-I-f + f + i)l (-! + 1)2+/. (-« - C + A;)2+a+h- 



(-4+a-fc)/2 

E ( 

h=0 



-ir 



3 I a _|_ c 

'2 "•" 2 "•" 2 



/i-^)l+;,(l)2+/.-i 

Analogously, the second and fourth summand of ( [4. 23] ) give 



(-4+a-fc)/2 



E 

h=0 



(-1) 



( 1 + - 

a+h-j V ^ 2 



' 1 (-§ + 1)2+/. i-a-c + k)3+a+h^ 



3 I g I e 

"I" 2 "•" 2 



(4.24) 



(4.25) 



Again we want to change the lower bound of the remaining single sums in ( 4.24[ ) and ( [4.25 ) 
to j — 2, respectively to j—3, according to the factor l/{l)2+h-j, respectively l/(l)3+/j_j. This 
is because then the two sums are 2-^1 -series and this makes it possible to apply Vandermonde's 
summation formula ( [4.22| ) again. In case of ( 4.25|) an d j = 2 this change is problematic since 
J — 3 < for j = 2. But the right-hand side of ( |4.16|) compensates this missing term. 

In terms of hypergeometric notation, it remains to show that 



(-ir(^ + l).(-a 



k)a 



V2 ^ 2 ^ 2 



k)il + ^ + l-j-kU+j 



2F1 






2 2 

+ 2F1 



§ + l+j) 



-1-f+j + i 



-c + k 



.S> T^ 2 ' 2 J 



k) 



0. 



2 f f + J + ^ 

For j > [a — k -\- l)/2 this is obvious because of the factor ((—a -|- k)/2)j. Otherwise, we 
use again Vandermonde's summation formula (4.22) with n = (a — 2j — k)/2, respectively 
n = (2 + a - 2j - k)/2. 

Finally we have completely proved that (6 -|- A;)^ is a factor of det(Da(6, c)). 

g 

re 5., 4. — n(& + c + l)i-i/((6 + c + 2)/2)p(a_2)/2l is a factor of det{Daib,c)): In order 

i=2 

to prove that the third and the fourth type of factors (see the beginning of Step 2) divide 
det(Da(6, c)), we could use the same procedure as for the factors of form (c + k) and the 
factors of form (6 + /c)^. But, fortunately, we will see that the triple sum that is the subject of 
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Section m provides an easier proof. 'Fortunately,' since in this case the procedure leads to quite 
complicated hypergeometric identities. Furthermore, the underlying linear combinations are 
more difficult to guess because an additional parameter originating in the higher multiplicity 
of the factors is involved. 

If we look at ([0| ) and (|45|), we see that ( |4^ ) is equal to ( |0| ) evaluated at x = (a + 6)/2 
and y = (a + c — l)/2 modulo some factors. Therefore a combination of ( ^.3D , ( |4.5D and 
Lemma |2| gives 



det{Da{b,c))= [f[(6 + c+l)i_i(i-l)! 

\i=2 
a a m 
X 

n=l m=l s=l 




b-a + 2 



EEE 



(-1) 



n+s 



l-a+c+2n \ 

2 /(a+fe-2)/2 



a-1 
a+b+2s 



a-1 



3— a+c '\ 

2 ){a+b-2)/2 V 2 ;(a+c-l)/2 



/ 2-a+fe '\ 
I 2 Jla 



(a+c-l)/2 (m-1 
-1 



(6 + c+l)s-i (n-1)! (m-n)! 



(4.26) 



We are aiming to show that nr=2(^ + c + l),j_i/((6 + c + 2)/2) [(a_2)/2] is a factor of (|4 
Using ( [4.26 ), it remains to show that the following sum is equal to a polynomial in b and c 
for all integers b and c if we fix a: 



^i=2 




6 + C + 2 



a + 3 



EEE 

n=l m=l s=l 



-1) 



r(a-2)/2] \ ^ / a-1 

/ l-a+c+2n N 
n+s^ 2 ;(a+6-2)/2 



6-a + 2 
2 

-a±b±2s\ 



a-1 



'3-a+c^ 



' (a+fe-2)/2 

{b + l)s-i 



(a+c-l)/2 / m. 
S - 



' 2-a+b \ 

. 2 /'(a+c-l)/2 

(c+ l)„_i (6 + c + n), 



(6 + C+l) 



n 



1)! 



(rn, — n)! 



(4.27) 



As it stands, this does not appear as a polynomial in b and c, since b and c occur in the 
second arguments of some Pochhammer symbols. But by using 



{s + n)r 

(s)r-j 



(4.28) 



which is valid for all integers n, we are able to ban b and c from the second arguments of the 
Pochhammer symbols in ( [4.271 ). Then, after some cancellations, we obtain 



a 



a m 



nc'-D'EEEM 



n+s 



i=2 



n=l m=l s=l 
' m — 1 



1- a + c + 2n 



b + c + l 



n-l 



-a + b + 2s 



X 



1 



'b±c±2\ ( b+c+l \ ,, ^, / -,^ /, n 

^ 2 )l{a-2)/2']\ 2 )s-l{b+l)s-l[C+l)n-l[b + C + n) 



' + C+l)s-l 



{n — l)!(m — n)\ 



for ( [4.27| ), and this is manifestly a polynomial in b and c. Indeed, because of s < m, < a, 
thanks to the binomial coefficient, the term (6 + c+ l)s-i cancels with the two terms on top 
of the fraction. 

Therefore we have finally proved that the linear factors in ( [4.9| ), respectively in ( 4.1[1| ), 
divide the determinant det{Da{b, c)) as a polynomial in b and c, and we may now turn to the 
computation of the remaining irreducible polynomial. 
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Step 3: Computation of the irreducible polynomial. We emphasize once more that 
Step 1 and Step 2 show that 



det{Da{b,c)) = (Linear factors in (|4.9| ), respectively ( 4.1C| )) x Pa{b,c), 



(4.29) 



where Pa{b, c) is a certain polynomial in b and c if we fix a. In this final step we prove that 
Pa{b,c) equals ( 4.11 ), respectively ( 4.121) . 

The computation of the irreducible polynomial, which will be denoted by Pa{b,c), is done 
in the following way: We consider the polynomial Pa{b, c) as a polynomial in b over Z[c] and 
find that it has 'nice' evaluations at b = —c — k, k odd and 1 < k < a, see ( 4.3C1| ) and ( [4.31 ). 
We are able to prove these 'nice' evaluations with the help of Lemma ^. Furthermore, we will 
see that Paib, c) is a polynomial in b with a degree smaller or equal to [(a — l)/2j . Therefore 
these 'nice' evaluations are just enough so that we can compute Pa{b, c) by using Lagrange's 
interpolation formula, see ( [4.32 ), respectively ( 4.33| ). 



First we will convince ourselves that the assertion about the degree of Pa{b, c) in b is true: 
The entry in the ith row and the jth column of the matrix Da{b, c), see ([4. 61 ), is a polynomial 
in b of degree j + 1, since H is a polynomial in b of degree 2. Hence, in the defining expansion 
of the determinant, each summand has degree 2 + 3 + --- + a = (a + 2)(a — l)/2 in b and 
therefore the degree of the determinant itself is at most (a + 2)(a — l)/2. One easily checks 
that the product of the linear factors is a polynomial in b with degree a^/2 in case of a is 
even and (a + l)(a - l)/2 otherwise. Because of that (a + 2)(a - l)/2 - a^/2 = [{a - l)/2j 
respectively (a + 2)(a — l)/2 — (a + l)(a — l)/2 = [(a — l)/2j is an upper estimation for the 
degree of Pa{b,c) in b. Later we will see that this is in fact the exact degree. 



We claim that 



Pa{-c-k,c) 



a-l 



n(^-i)! 



\i=2 




(^)!(_l){^-i)/2 

f a-k+l \ 
(a-l)/2 I, 2 J{k-l)/2 



c + k+1 



(a-fc)/2 



C+1 



(4.30) 



(fe-l)/2 



if a is odd and /c = 1, 3, . . . , a, and 



Pa(-c- k,c) 




(Y)!(-l)(^-^)/^ 

{a-2)/2 l^-j(A:-l)/2 

c + k + 2\ fc + 2 



(4.31) 



(a-fe-l)/2 \ ^ / (fc-l)/2 

Assuming the truth of the claim, we would have 

c+l\ /l + b + c\ /c + 2k + 2^ 



if a is even and k = 1,3, . . . ,a — 1. 



a-l a 



Pa{b,c) = (-] Uii-iy. J2 

^ ^ i=2 A:=0 



6 + C + 2/C + 3 



k \ ^ / (a-2fc-l)/2 

(2)(a-2fc-l)/2 



(a-2fc-l)/2 (l)(a-2fc-l)/2 



(4.32) 
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if a is odd, and, respectively, 



(a-2)/2 



^a(6,c)= - n(^-l)! Y. 

^ ' 1=2 k=0 



c + 2\ fl + b + c\ fc + 2k + 3' 

k \ ^ / k \ ^ J {a-2fe-2)/2 

b + C + 2k + 3\ (2)(a-2Jt-2)/2 



(a-2A:-2)/2 (l)(a-2fe-2)/2 



(4.33) 



if a is even, by Lagrange interpolation. 

Thus, it remains to show (|4.30| ) and ( 4.31 ). By definition, Pa{b,c) is the quotient of 



det(-Da(^5 c)) and the linear factors in (^^), respectively in ( 4.10[) (see ( [4.29| )). This, together 
with a tedious but straightforward calculation shows that the claims in ( [4.301 ) and ( [4.31| ) 
conflate to the following 



det(Z),(6,c))( J_, 



b+c-l 

2 
g+c— 1 

2 



[](l + 6 + c),_i(z-l)! 



\i=2 



a-1 



b-a + 2 



a-1 



(-1) 



a(c+i; 



fc-1 



-c—k 



(fc-1)! 

(4.34) 



The reader should notice that we are not able to directly set b = — c — A; on the left-hand side 
of ( |4.34 ), because the denominator vanishes for b = —c — k. As mentioned before we want to 



replace the determinant on the left-hand sides of (4.34) by the triple sum in Lemma g. Still 
even after cancellations we are not able to directly set b = —c — k in this triple sum, since 
then the denominator of the summand becomes zero if A; > n. 

In order to avoid an indefinite expression we act as follows: If we reexamine the proof of 
Lemma |2|, we obtain the following generalization of Lemma |2|: 

Let a,b,c he positive integers and b' ,x,x he integers. Then the determinant 



det 

l<«.i<a+l 



/ 


1 < J < a 


j = a + l 

/b'-x+y\ 
\y-i+l) 


\ 


r7-V) 






1 < j < a 



a + l 



(4.35) 



equals 

a b c 

nnn 



i + j + k-i\ (1), 



^1=1 j=l k=l I 

'" ' c^ X — y ^n — 2\ fb' — a;-|-y-|-s — 1^ (6+1 

X — 1 /\ 6' — x + s — 1 



X 



EEE 

n=l m=l s=l 



(-1) 



n+s 



ls-1 



6 + c + l),_i 
m - 1\ (c + l)„_i (6 + c + n)m- 



s — lj (n — 1)! {m — n)\ 



. (4.36) 



We set b = —c—k + e and b' = —c—k in (4.36). In view of the fact that ( [4. 351) is continuous 
in 6, and if we compare the determinant in (4.6) with the determinant in (p. 6]) evaluated at 
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X = (a + 6)/2 and y = (a + c — l)/2, we obtain that the left-hand side of ([4.34|) is equal to 



a a m 



71=1 m=l s=l 



-1) 



n+s 



-fc+2n-3 \ / -k+2s-3 
2 



a—c—k—2 



2 

g+c— 1 
2 



m — 1 
s- 1 



(-C - A: + g + l)^_i (c + l)n_i (-fc + e + n)r 
{-k + e + l)s-i (n-1)! {m - n)l 



(4.37) 



Now it remains to show that the right-hand side of ( [4.34 ) equals ( 4.37| ) for k = 1,3, . . . ,2[a — 
1/2J. 

Therefore let us consider the triple sum ( [4.37| ). First of all, we are allowed to extend the 
sum over s to all positive integers, since {"^Zi) = if s > m. Using hypergeometric notation 
for this innermost sum, we get 



lim > > 

71=1 711=1 



(-1) 



n+1 



3-^2 



1 — m, 1 — c - 
l-k + e,l 



K -\- E, 2 



k 

2;1 



'1 + c)_i+„ (1 - § - f - |)_ i+a+£ (I - 5 + I + "')-i4-£_£-*L {-k + e + n 



^l)-i + £. + £ (l)_i-L£_£_fe (l)-m-7i (1)-1+^ 
2^2^2 ^2 2 2 



We apply the following transformation formula due to Thomae ||5|, (3.1.1)] to this hypergeo- 
metric sum. 



3-^2 



a, 6, — n 



(e)n 



' 7? 
-3-C^2 



—n,b, —a + d 
d,l + b — e — n 



;i 



(4.38) 



with a = l-m, 6 = l-c-fc + e, n=(/c- l)/2, d = 1 - fc + e and e = 1 - a/2 - c/2 - k/2. 
Writing the resulting 3i^2~series as a sum over s, after some simplifications and cancellations 
this gives 



lim 

e 



a a oo 

jEEE 

71=1 771=1 S=0 



(-1) 



i-^ + f+fc + 71+S 



a c k 

-2 + 2 + 2 

-1 a 



1 + c)_i+„ (1 - c - A: + e)s (^ - f + f + n,)_-^^a_c_fc (-A; + e + 



n 



(l)m-7i (l)-l+n (l)_l I fc_s (l)s (l)l+S: 



2^"2 2 "'^^ 



(l-A: + e) 



Next we interchange the two inner sums and reverse the order of summation in the inner- 
most sum. We obtain 



a oo — 1+a 

l-isEEE 

71=1 s=0 m=0 



(-1)5 



]+k+n+s 



a c k 

'2 + 2 + 2 



(1 + c)_i+„ (1 - c - A: + e), (i - f + § + n)_i^._c_| i-k + e + n)^_^_^_^^ 

il)a-m-n (l)-l+n (l)_i + fc ^ (l)s (1) i + "_§_fc+. (1 " A: + e)s 



The reason why reversing was so helpful, is that we now can extend the summation with 
respect to m to all nonnegative integers (since l/{l)a-m~n = if m > a — n). If we use 
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hypergeometric notation this yields 



lim > > 



e-^0 



n=l s=0 



(-1)^ 



}+k+n+s 



2F1 



1, —a + n 

a + k + s ■ 



(1 + c)„i+„, (1 - c - A; + e), (^ + § + I) (i - f + f + n)_,^„_c fe (-/. + £ + n)„_„, 



(!)„_„ (l)_i+„(l)_i,._fl),(l). 



_^^^{l-k + e). 



We apply Vandernionde's summation formula (4.22) to this hypergeometric sum. After 
some manipulations this gives 



J^nEE 



£-♦0 



n=l s=0 



-1) 



a 
'2 



c k 

2 + 2 



-a + k - e - s)a-n {1 + k - e 
(l)i+|_„(l)_i+^_. 



n 



(l + c)„i+„(l 



C fC + £)s (22 2 



fc + ■s)-l+A;-n-s 



(l)a- 



1 



-1+n 



(l),(l-fe + e) 



we 



Now we are able to perform the limit e — > 0, since, because of the factor l/(l)-i/2+fc/2 
have k > s. 'Performing the limit' means that we simply set e = 0. 

With pleasure we notice that our triple sum has simplified to a double sum. Big surprise 
arises when I finally claim that the summand of the sum is only different from zero if and 
only if n = k/2 + 1/2 and s = k/2 — 1/2, and we therefore get rid of all sums: The factor 
l/(l)i/2+A;/2-n implies n < 1/2 + k/2, and the factor l/(l)_i/2+fc/2-s implies s < —l/2 + k/2. 
Next we notice that 



n 



(k — n — s) 

The numerator on the right-hand side is zero since —a + k — s < —a + k < and k — n >0. 
Therefore the summand is only different from zero if A: — n — s = 0. This together with our 
first observation gives my assertion. 

It remains to compute the summand in the last double sum for n = k/2 + 1/2 and s = 
k/2 — 1/2 (and, of course, e = 0). After some simplifications one does indeed get the right- 
hand side of ( [4.34| ) and Lemma |^ is finally proved. 

Now we turn to the case that the side lengths a,b,c of the hexagon have the same parity 
and therefore to the evaluation of the determinant in Lemma || with x = {a + b)/2 and 
y = (a + c)/2 (see the beginning of this section). 

Lemma 4. Let a,b,c he integers and a = c (mod 2). Then the determinant 



det 



/ 


1 < i < a 

( ^ + ^ ^ 
Vc - i + j) 


j = a + l 


\ 


/ b+c-2 \ 

v§-4+J 






1 < i < a 



i = a + 1 



is equal to 
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b c 



nnn 



i+j+k-l 



(l)c 



lfc=V + J + ^-2 ) {b+l)c+a- 



c-2\ /a+b+c-l 



b-1 



'c+ 1\ /6 + C + 2 



b-1 
2 



2 \ Q''-! 



(a-l)/2 



+ E 



2/'(a-l)/2 
2 ''(a-l)/2 V 2 /(a-l)/2 (l)(a-l)/2 



.2J(a-l 



c + 2 



fc^i V 2 yfc_iV 2 /fc 



6 + c\ /c + 2A; + l 



(a-2fc+l)/2 
b + C+2k + 2\ (2)(a-2fc-l)/2 



(a-2fc-l)/2 (l)(a-2fc-l)/2 



in case that a is odd, and 







C+2 



b + c + 2 



^2)a/2 



2 ^(a-2)/2V 2 /a/2 (l)(a-2)/2 



c + l\ /6 + c\ /c + 2A; + 2 



fc=l \2/^, V2/^\ 2 / (a-2A;)/2 



6 + c + 2A; + 2 



(2)(a-2fc)/2 
(a-2A:)/2 (l)(a-2fc)/2 _ 



in case that a is even. 



As already mentioned at the beginning of this section the proofs of Lemma ^ and Lemma Q 
are quite similar. Therefore we restrict ourselves to just pointing out the differences to the 
proof of Lemma ^ but we omit the details. The only major difference to Lemma |3| is that 
in the present case the irreducible polynomial Pa{b,c) (see, e.g., ( [4.29| )) has one 'exceptional 
evaluation' which has to be treated in a different way than the rest. By the way, this 
exceptional evaluation also causes the formulas for the number of rhombus tilings which 
contain the 'almost central' rhombus above the centre to be not as simple compared to the 
formulas for the number of rhombus tilings which contain the central rhombus. 

Proof of Lemma ^ 

Step 1: Prom our determinant to a determinant with polynomial entries. 

We can straightforwardly adopt Step 1 from the proof of Lemma y, where we have reduced 
the problem to the evaluation a polynomial determinant, since we have not even specialised 
the coordinates of the fixed rhombus {x, y) there. In the following Da{b, c) denotes the matrix 
underlying the determinant in ( ^^ ) evaluated for x = {a + b)/2 and y = {a + c)/2. 

Again I was led to a conjecture concerning the determinant det(Z)a(6, c)) by computing it 
for small values of a: 
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^ ' (a-l)/2 V 2 / (a-l)/2 V ^ y (a-l)/2 

,,C+1^^ /^6 + C + 2^^ (2)(a-l)/2 

X 



2 y{a-l)/2 V 2 /(a-l)/2 (l)(a-l)/2 

j^^l \ -^ / fc-l V -^ / A: \ ^ / (a--2A:+l)/2 



6 + c + 2A; + 2^ (^)(a-2fc-i)/2 

{a-2fc-l)/2 (l)(a-2A:-l)/2^ 



if a is odd, and 



- \i=2 / \i=2 

b\ fh + 2\ fc + 2\ f2 + b + c\ ^a-i 



^'a/2\ 2 / {a-2)/2 V ^ / (a-2)/2 V ^ / (a-2)/2 

' 'c + 2\ /6 + C + 2\ (^)a/2 



2 /(a-2)/2V 2 /a/2 (l)(a-2)/2 

f^^lX-^/kX-^/kX ^ J {a-2k)/2 



h + C + 2k + 2\ {\){a-2k)/2 

{a-2k)/2 (l)(a-2fc)/2^ 



(4.39) 



(4.40) 



if a is even. Just as in the proof of Lemma |3| the two assertions ( 4.39] ) and ( [4.40| ) are equivalent 
to the assertion in Lemma ^. 

The hnear factors of det(L'a(6, c)) are 

n(l,,,eu(^)^ (£±i) (^±|±£) 2-^ (4.41, 

j=2 V -^ / (a-l)/2 V ^ / (a-l)/2 V ^ / (a-l)/2 

if a is odd, and, respectively, 

j=2 V^/a/2V ^ / (a-2)/2 V ^ / (a-2)/2 V ^ / (a-2)/2 
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if a is even. The irreducible polynomial, denoted by Pa{b,c), reads 

ly-^ f fc+l\ fb + C+2\ (l)(a-l)/2 



a 
\i=2 



m 




(a-l)/2 V 2 ) (a-\)j2 (l)(a-l)/2 



fe-l \ ^ / A; \ ^ / (a-2fc+l)/2 

6 + c + 2A: + 2^^ {\)(^a-2k-\)l2 

{a-2k-l)/2 (l)(a-2fc-l)/2^ 



(4.43) 



if a is odd, and respectively, 

'c + 2' 






/ 

-2)/2 V 


6 + c+2\ (i), 


1/2 






(a- 


2 A/2(l)(a- 


2)/2 




\ 


(h + 


c^^ l^c + 2A; + 2^^ 








/fc 


V 2 


/ fc V 2 y („ 


-2fc)/2 








^^ ^6 + c + 2A; + 2^^ 




(|)(a- 


-2fc)/2 






V 2 ;, 


i-2fc)/2 


(l)(a- 


-2A:)/2 



(4.44) 

V ^ / (a-2fc)/2 y''){a-2k)/2 J 

if a is even. 

Step 2: The linear factors in ( 4.41| ), respectively in ( 4.42| ), divide det(Da(6, c)) as 
a polynomial in h and c. 

Again we have four different types of linear factors: 

1. Factors of the form [c + k): If a is odd, k = 1,3, ... ,o — 2. Otherwise A: = 2,4, . . . ,a — 2. 

2. Factors of the form (h+kY: If a is odd, /c = 1,3, . . . ,a — 2. Otherwise A; = 2,4, . . . ,a — 2. 
The factor h occurs once if a even. 

3. Factors of the form (6 + c + kf-'^-^ for /;: = 1, 3, . . . , 2[(a - 2)/2] - 1. 

4. Factors of the form (6 + c + kf-'' for A; = 2, 4, . . . , 2 [(a - 2)/2j . 

re 1. — The factors of the form {c + k) divide det{Da{b, c)): As described in the analogous 
passage of the proof of Lemma y, each linear factor of det(Z?a(6, c)) corresponds to a linear 
combination of rows of a certain matrix. The linear combinations for the factors of the form 
(c + k) are 

a-fc+l (_lY-i(U<j) ■ ( -a+k-2+2i \ 

E ' ' ,^ il±-.; '°--^- P.(^-^(..,=0 (4.45) 

i=(a-k+2)/2 \^)a-k+l-i { 2 J a-k+l-i 

for j = 2, 3, . . . , a, 1 < k < a — 1 and k = a (mod 2). In order to prove that (c + k) divides 
the determinant det{Da{b, c)) we just have to prove the identity ( 4.45 ). 

The comparison of the identity in ( 4.13 ) and the identity in ( 4.45 ) shows that these two 
identities are quite similar. And in fact the proofs do not differ essentially from each other, 
either: Again we are able to apply Gosper's algorithm ^ for hyper geometric sums to the 
left-hand side of ( [4.45 ) and recognize that this is actually a telescoping sum. 

The expression analogous to ( [4.14| ) is 

^ (-l)J' (1 - f + I) . (-a - 6 + fe) .+^._2 

2^ 77^; r^ a b~T\ ^\x=(a+b}/2,y=(a+c)/2,c=-k,i^a-k+l-i " 0' (4-46) 
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which is the identity in ( [4.45 ) after reversing the summation order. If f(i) denotes the 
summand in the previous sum and g{i) denotes the following expression 



(-l)^(l-f + |)J-a-6 + A:) 



i+j-2 



(1)- 



l-a+i+i 



+ ^).-l 



h-hi + 2j + bj -bk + ik- jk + k"' 



then we have 



f{i)=g{i + l)-g{t). 



By using this identity it is easy to compute the left-hand side of ( 4.46 ) and to show that it 
is equal to zero. 

re 2. — The factors of the form (b + k) divide det{Da{b, c)): 

The two linearly independent linear combination for the factors of the form are (b + k) are 



{a+k+2)/2 

E ■ 

i=k+2 



ly-^-^ic + a-i + l) 



i—k 



'a+fc-2^+4^ 



~k-2 



(1) 



i-k-l 



'c+a-2i+4^2 



i-k-2 



Pi^k-i{c + a-k-i + 2)Da{-k,c)^ij^ =0 (4.47) 



ifO<A;<a — 2 and k = a (mod 2) , and 



(''+!^y^{-iy-''{c + a-i + 2 

j=fc+3 



'a+fc-2^+4^ 



(1) 



c+a-2i+4A2 



J-fc-1 V 2 n-k-2 

X Pi-k-2{c + a-k-i + l)ba{-k, c)(^ij) = -Da{-k, c)(fc+ij) (4.48) 



ifO<A;<a — 2 and k = a (mod 2), where Pnic) is the same sequence of polynomials as in 
Lemma ^, see (|4.17 ). The first identity, ( 4.47 ), corresponds to ( 4.15|) , and the second identity, 
( [4.48 ), corresponds to ( [4.16| ). Again these two linear combinations do not cover the case that 
k = a — 2. But similar to the situation in Lemma y it can be verified directly that (6 + a — 2)^ 
is a factor of the determinant det{Da{b, c)). 

The proofs of the identities ( 4.47 ) and ( [4.48| ) are analogous to the proofs of their cor- 
responding identity: Again we first have to interchange the summation and then split the 
double sum into four smaller sums according to the polynomial H. Just as in Lemma |3| the 
inner sums of the first and the third double sum nearly cancel out each other, and so do the 
inner sums of the second and the fourth summand. This reduces the problem to identities 



only involving 2-Fi-series and therefore Vandermonde's summation formula (4.22) finishes the 
proof of these identities. 

a 

re 3., 4- — Y\{b + c+l)i-i/{{b + c+l)/2)^(^a-i)/2'] is a factor of det{Da{b,c)): Just as in 

i=2 

Lemma y the triple sum from Lemma g provides an easy proof that the factors of type 3 and 
type 4 divide the determinant det{Da(b, c)) as a polynomial in b and c. Again a combination 
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-a + 2\ fb-a + 2'^ 
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d&i{ba{h,c))= rQ(6 + c+l)i_i(i 

\i=2 
a a m 
X 

n=l 971=1 s = l 




EEE 



a+c+2ra^ 



-1 



in+s 



'(a+6-2)/2 



a-1 
a+6+2s^ 



a-1 



V 2 /(a+fe-2)/2 



'2-a+6^ 



'(a+c)/2 /^rn 
S - 



(6 + i; 



s-l 



V 2 7(a+c)/2 

(c+l)„_i(6 + c + n) 



f4 491 
{h + c+l)s-i (n-1)! (m-n)! J' ^- ^ 

By using ( [4.28| ) we obtain that the quotient of ( |4.49D and the product of the linear factors 

a 

]^ (6 + c + l)j_i/((6 + c + l)/2) r(a-i)/2l is the following polynomial in h and c if we fix a: 



1=2 



m-^y-EEEi-^r- 



i=2 



n=l m=l s=l 



-a + c + 2n 



b + c 



n-l 



-a + b + 2s 



m 

s - 



1 



hc+l' 



' i{a-i)m 



'h+c+2 



)s-i {b + l)s-i(c + l)„_i(6 + c + n)r 



ly {b + c+l)s-i {n-l)\{m-n)\ 

Step 3: Computation of the irreducible polynomial. Analogously to the situation 
in Lemma ^ we now evaluate the quotient of the linear factors in ( 4.41 ), respectively ( 4.42 ), 
and the determinant det(Da(6, c)). This quotient, which will be denoted by Pa{b,c)^ is a 
polynomial in b and c if we fix a. Again we will find that the polynomial Pa{b, c) has enough 
'nice' evaluations so that we can compute the polynomial by using Lagrange's interpolation 
formula. Namely, these 'nice' evaluations arise for b = —c — k, < k < a and k is even, and 
the degree of Pa{b,c) as a polynomial in b happens to be [a/2j. 

The assertion about the degree of Pa(b,c) as a polynomial in b can be checked routinely, 
just as in Lemma & Concerning the computation of 'nice' evaluations of the polynomial 
Pa{b, c) there arises the one and only major difference to Lemma ^ For k = 2,4,..., [a/2j 
the evaluations Pa{—c — k,c) can be proved in the same way as in Lemma y, but for A; = 
the situation is different, and therefore this exceptional evaluation needs a separate proof. 
This fact already shows up if we look at the following conjecture for the 'nice' evaluations: 

We claim that 



Pai-C,c) 

if a is odd, and 

Pai-C,c) 



a-1 



\i=2 




(a-l)/2 



C+1 



(a-l)/2 



^)"ln.-.-)(f)G)J^) 



(a-2)/2 



(4.50) 



(4.51) 



if a is even, and, furthermore, 

a-l / I 



Pa{-C- k,c) 



m 

\i=2 




;|)!(-l)^/2 



(a-l)/2 

c+k + 1 



a—k \ 
2 )k/2 



{a-k+l)/2 



C + 2 



(fc-2)/2 



(4.52) 
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if a is odd and A; = 2, 4, . . . , a — 1, and 



Pa{-C-k,c) 



-2 a-1 

n 

i=l 



(|)!(-l)fe/^ 

( a-k+l\ 
(a-2)/2 \-^—){k-2)/2 



c + k + 2 



{a-k)/2 



C+1 



fc/2 



(4.53) 



if a is even and A; = 2, 4, . . . , a. 

Assuming the truth of the claim, we would have 



Pa{b,c) 



a-1 / a 

11(^-1)' 
\i=2 

(a-l)/2 

+ E 

fc=i 



c+1 



b + c + 2 



(a-l)/2 

c + 2\ /6 + c\ /c + 2A; + l 

fc-i 



(2)(a-l 



{a-l)/2 



(a-l)/2 (l)(a-l)/2 



(a-2fc+l)/2 



X 



b + c + 2k + 2 



(2)(a-2fc-l)/2 
(a-2fc-l)/2 (l)(a-2fc-l)/2 



if a is odd, and, respectively 



Paib,c) 



1 X a-2 a I / , r, 



4=2 



6 + C + 2 



^2Ja/2 



a/2 

+E 

fc=i 



(a-2)/2 

c + l\ fb + c\ fc + 2k + 2 



a/2 (l)(a-2)/2 



{a-2fe)/2 

6 + c + 2fe + 2^ 



(2)(a-2fc)/2 

(a-2fc)/2 (l)(a-2fe)/2 



if a is even, by Lagrange interpolation. 

Thus, it remains to show (4.50), (4.51), ( ^.53| ) and ( |4.52| ). In these claims we replace 
the polynomial Pa{b,c) by the quotient of det{IDa{b,c)) and the linear factors in ([4.41| ), 
respectively ( 4.42| ), in order to see that the four claims conflate to the following two: 



• h+c-2 \ / b+c 

2det(D.(6,c))(^J_,j(^4^ 



J{{l + b + c),.i{i-l)\ 



c-a + 2 



a-1 



b-a + 2 



(b + c) 



a-1 



b=-c 
a+1 I I / i\a+l 



-a + (-l)°+^a + c + (-ir+^c 
c {a + c) 



, (4.54) 
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and 



b+c-2 



det{D,{b,c))(J^ 

2 ^ ^ 2 



b+c 
2 



f];(l + 6 + c),_i(i-l)! 



\i=2 



c- a + 2 



a-\ 



b- a + 2 



a-\ 



(-1) 



a-1 



(c + r 



fc-i 



{k-\)\ 



-c—k 



(4.55) 



for k = 2,4,..., [a/2j . Again the reader should know that we are not able to directly set 
6 = — c — A; on the left hand sides of ( 4.55 ) and ( |4.54 ), because the denominator vanishes for 
b = —c — k. 

First we consider the case that k ^ 0. As already mentioned, the situation in this case is 
quite the same as in Lemma y. T herefore I only explain the essential steps and omit details. 

For the left-hand side of ( [4.55| ) we use again the modification of Lemma |2|, see ( 4.35 ) and 
( 4.36| ). Thus, this left-hand side is equal to 



n=l m=l s=l 



-1) 



n+s 



-fc+2ra-4 



a— c—k— 2 



-k+2s-2 

2 

a+c 

2 



m — 1 
s-1 



(-C - k + e + l)s-i (c + l)n-i {-k + e + n)r 



-k + e + l)s-i 



(n-l)l 



(m — n)! 



. (4.56) 



Hence, it remains to show that the right hand side of ( 4.55| ) is equal to ( 4.56| ) for k = 
2,4,...,La/2j. 

Just as in Lemma |2| we apply the transformation formula due to Thomae, see ( [4.38 ), to 
the innermost sum of the triple sum, with slightly changed parameters, i.e., a = 1 — m, 
b=l — c— k + e,n = {k — 2)/2, d = 1 — k + e and e = 1 — a/2 — c/2 — k/2, compared to the 
proof of Lemma 0. This yields 



a a oo 

jT„EEE(-i) 

n=l m=l s=0 



'. + i,+n+s 



-a c k 

Y^2^2 



(l + c)_i+„(l-c-A:),(^ + f + n) 



-1+1 



(e — k + n) 



(!)„_„ (l)_i+„ (l)_i+fc ^ (l)s (l)i+f-f-fe+s (1 + e - k)s 



Next we interchange the two inner sums, reverse the order of the summation in the new 
innermost sum and use hyper geometric notation for this innermost sum. We obtain 



lim 

e 



m oo 



ra=l s=0 

(l + c)_i+„(l-c-A;) 
X 



1, —a + n 

1 — a — e + k 



« V 2 ' 2 ' 



2M 2" + 2 +"')-l+§-f 2 



k{e-k + n)a-n+s 



(!),_„ (!)_!+„ (l)_,^._^(l),(l)l+. 



-k+s 



(l + e-k). 
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Then we apply Vandermonde's suniniation formula ( 4.22| ) to this hypergeoinetric sum. After 
some manipulations this gives 



SEE (-')-' 



e->0 



n=l s=0 



Y + l + ^l (l+c)-l+n(l-C-/c), 



-a- e + k- s)a~.n (1 - e + /c - n - s)s (2 + f 



k + s)-i+k-n^ 



(l)a-n (l)i+|_„ (l)-l+n (l)_i+|_, (1). (l + S - k)s 

Now we are able to perform the limit e — > and then observe that the summand of the 
double sum is only different from zero if n = k/2 + 1 and s = k/2 — 1. The evaluation of the 
summand of the double sum for these special values of n and s finishes the proof of ( [4.55| ) 
for A; = 2, 4, . . . [a/2j . 

Finally we consider the case that k = 0, see ( [4.54| ). By Lemma |2|, ( [4. 3D and (^]^) the 
left-hand side of ( 4.54 ) is equal to 

2 ^^ -(_!)"+- (i + 6)_i_,,(l + c)_i+„ 



b + c 



71=1 m=l s=l 



(1)^ 



1+s 



-1+n 



;i + f + §)_i_<i + i>+^ (^ + § + ^)_l+ii + |i {b + c + n)m-n {l + m 



S)~~l+s 



(l)_l+a+6 (l)m-n (l)_i_£+6+s {l+b + c)_i+s 

before evaluating it at 6 = — c. We take the factor (6 + c)/2 out of the Pochhammer symbol 



2 + a + c 



(-2-a+b+2s)/2 



2 + a + c 



(-2+a+b)/2 



6 c 

2 "^2 



6 + C + 2 



and set b = —c. After some cancellations and simplifications concerning the summand of this 
triple sum, this yields 



a a m 

EEE 

n=l m=l s=l 



(-1)" 



-l+t~g+n+s 



(1 -c)„i+s(l +c)_i+„ 



J)-l+n (l)-l+s 

^ (f - f )l-n (1 + f + i)-l-f-§ Nm-n (l + m - s)_i+, 
(l)l-n (l)m-n (l)-l-S_£+s 



(4.57) 



The factor l/(l)i_„ shows that the summand of the triple sum is only different from zero if 
n = 1. Therefore we get rid of the sum over n. We use hyper geometric notation for the inner 
sum and obtain 



E 

m=l 



1 — c,l — m_ 

1 a c ) 

,2 2 



a I c 
2 "•" 2 



Next we apply Vandermonde's summation formula ( 4.22| ) to the 2-^i^series. This yields 



E 

m=l 



V 2 llf 



-l+m 



(4.58) 
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Once more, Gosper's algorithm Q can be used to see that this single sum is actually tele- 
scoping. Namely, there holds 



.Nc V 2 + 2)-!- 
^ (-1- 



a c 

'2 2 



+ m){-l-l + 



2/m+l 



-f-f+m-l)(-l-f + f), 



U + 2 + 2/ "^ V I 2 2ym+l V 2 2/m / 

Because of that, and because a and c have the same parity, the single sum ( |4.58 ) is equal to 

-a + (-l)"+^a + c+(-l)°+^c 
c {a + c) 



If we compare this to the right-hand side of ( [4 .541 ) we see that this is just what we claimed. 
This finally completes the proof of Lemma 0. 



5. The Proofs of Theorem || and Theorem || 

We end this article with the proofs of Theorem |3| and Theorem Q. We start with Theorem ^. 
In Theorem ffl we showed that the probability to choose a rhombus tilings of a hexagon 
with side lengths a,b,c,a,b,c which contains a rhombus in the centre is, using hypergeometric 
notation. 



4F3 



1 Q:i_l_£i_l_^_l_£l 

2 2' 2 "•" 2' 2 """ 2 "•" 2' -*-. 1 
1 ^1_|_c 3 I b I c i-"- 



2' 



2' 2 + 2 + 2 



(^)_i+g(l)c (! + §)_ !,.(! + 



2^-l+2 + £ (2 + 2 + 2) 



X 



— 4- — 
" 2"'"2 



_1 I Q_ 



2a-l 



(l)-i- 



i+6 (1 -I- 6)_i+a_|.c 

2~2 



(5.1) 



in case that a is odd, and 
4-P3 



1_« 1_lC lifcic -|^ 



2'^ ' 2'2 ' 2 ' 2'.1 

3 a 3_i_c 3_i_fe_i_c' 

2 2'2"'"2'2"'"2"'"2 



X 



(i)_i+. (1). (6) (i + f ) . (I + f )_i+|+| (i + I + §)-i+|2 



a-2 



(l)_l+Wl)6^(l-f6). 



(5.2) 



-l+a+c 



in case that a is even. In order to see this, divide (|1.2| ), respectively (O), by MacMahon's 
formula for the total number of rhombus tilings of a hexagon with side lengths a,b,c,a,b,c 
given in ( |1.1| ). 

Next we apply Bailey's transformation formula (see p3| , (4.3.5.1)]) between two balanced 
4F3-series, 



4-f^3 



a, b, c, —n 
e,f,l + a + b + c — e 



1 



- a)n{f -a)r 

(e)n(/)n 



4-^3 



f -n 
-n, a,l + a + c — e — f — n,l + a + b - 



f -n 



1+a+b+c 



f — n,l + a 



n,l + a — f 



1 



n 



(5.3) 



which is valid if n is a positive integer, with a = 1,6 = (1 + c)/2,c = (1 -|- 6 -|- c)/2,n 
(a - l)/2, e = (2 - a)/2, f = {c + 2)/2 in case of a is odd, and with a = l,b = {2 + c)/2, c 
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(1 + 6 + c)/2, n = {a — 2)/2, e = (3 — a)/2, / = (c + 3)/2 in case of a is even. This gives 



4-^3 



1,1, 



3 3 I 6 I _ 

2 ' 2 "•" 2 """ 2 ' 2 



,1 + 

c 3 a 

2 



21 ' 2 -1 



a!6!c!((a + 6 + c - 2)/2)!((a + 6 + c)/2)! 



((a-l)/2)!2((6 
if a is odd, and 



l)/2)!2(c/2)!((c - 2)/2)!(a + b + c- 1)!((-1 + a + c)/2)((l + b + c)/2) 



4-f^3 



1,1,1 



3 3 I 6 _|_ c 3 

2 ' 2 "•" 2 """ 2 ' 2 



2' 2 "•" 2 



c;i 

2 



(a-l)!6!c!((a + 6- 



l)/2)!2 



((a - 2)/2)!2(6/2)!2((c - l)/2)!2(a + b + c- 1)!((-1 + a + c)/2)((l + b + c)/2) 

if a is even, after transforming the Pochhammer symbol (a)„ = {a + n — l)!/(a — 1)! to 
factorials. 

Now we substitute a ~ aN, b ~ f3N and c ~ 7A^ and perform the limit N -^ c«. We use 
Stirling's formula to determine the limit for the quotient of the factorials in the second line 
as ^/a\[]5^\Ja + /? + '^/{2T:{a + 7)(/3 + 7))- For the 4F3-series, we may exchange limit and 
summation by uniform convergence: 

1 1 i_a 1_L^ 



lim 4F3 



1 

2 2- 
3 3 I 6 I c 3 



2 ' 2 + 2 + 2 ' 2 



2-^1 



1,1 

3 ; 

2 



a/3 

(/? + 7)(a + 7) 



A combination of these results and the use of the identity 

"1,1 



2F1 



3 ;^ 
2 



arcsm v z 



V^ 



establishes Theorem 



The proof of Theorem ^ is analogous: Again we divide the formulas in Theorem g, see ( |1.4| ) 
and (|1.3| ), by MacMahon's formula for the total number of rhombus tilings of a hexagon with 
side lengths a, 6, c, a, 6, c given in ( |1.1| ), and obtain the probability to choose a rhombus tiling 
of a hexagon with side lengths a, 6, c, a, 6, c which contain the 'almost central' rhombus above 
the centre. Next we apply Bailey's transformation formula ( ^.31) on the hypergeometric sum 
and transform the Pochhammer symbols to factorials. We then substitute a ~ aN , b ~ PN 
and c ~ 7A^ and perform the limit A^ — > 00 in the same way as before. The additional 
summands (caused by the 'exceptional' evaluation of the irreducible polynomial Pa{b,c)) 
that appear in the formulas in Theorem y (when compared to the formulas in Theorem ||) 
vanish when this limit is performed. 
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A hexagon with side lengths a, b, c, a, b, c, where a = 3, & = 5, c = 4. 

Figure 1. 




A rhombus tihng of a hexagon with side lengths a, b, c, a, b, c. 

Figure 2. 
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A hexagon with side lengths a, b, c, a, b, c, where the central rhombus is marked. 

Figure 3. 




A hexagon with side lengths a, b, c, a, 5, c, 

where the 'almost central' rhombus above the center is marked. 

Figure 4. 
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A rhombus tiling of a hexagon with an indication 

of the corresponding family of non-intersecting lattice paths. 

Figure 5. 



Av^ 



A^»- 
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{x- l,y) {x,y) 



E3 



— • 
E2 



El 
The family of non-intersecting lattice paths corresponding to the rhombus tiling in Figure H. 

Figure 6. 
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The oblique angled coordinate system. 
Figure 7. 



